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ABSTRACT 


y  This  paper  considers  a  joint  replenishment  inventory  problem  with  a  continuous- 
review  (S,  c,  s)  policy  for  the  backorder  case  with  Poisson  demands  and  constant 
procurement  lead  times.  ^  *• 

Whenever  item  i's  inventory  level  hits  sj  (reorder  point)  or  lower  it  triggers  an 
order  so  as  to  raise  item  i's  level  to  S:  (order  up  point).  At  the  same  time  any  other 
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item  j  with  inventory  level  at-or- below  its  can- order  point  c-j  is1  included  in  the 
replenishment. 

A  Poisson  demand  model  with  a  queueing  description  of  the  system  s  operation  is 
analysed,  and  comparisions  are  conducted  for  joint  versus  individual  orders  in  the  case 
of  multi-item  problems,  where  joint  replenishment  of  several  items  may  reduce  setup 
costs.  <  ^  , 
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I.  INTRODUCTION 


A.  DESCRIPTION  OF  THE  PROBLEM 

The  purpose  in  constructing  a  mathematical  model  of  an  inventory  system  is  to 
use  it  as  an  aid  in  developing  a  suitable  operating  doctrine  for  the  system.  The 
criterion  most  frequently  used  for  selecting  the  operating  doctrine  is  that  of  profit 
maximization  or  cost  minimization.  In  some  cases  the  task  of  determining  the  optimal 
operating  doctrine  is  so  difficult  that  it  is  either  impossible  or  uneconomical  to 
determine  the  optimal  doctrine,  and  instead,  one  optimizes  with  respect  to  some  subset 
of  operating  doctrines.  Occasionally,  the  mathematical  model  may  be  so  complicated 
that  it  is  extremely  difficult  to  do  anything  analytically.  In  such  situations,  simulations 
are  used  to  study  various  operating  doctrines. 

We  are  concerned  with  finding  optimal  ordering  policies  for  a  multi-item 
inventory  system  in  which  the  inventory  position  is  under  continuous  review.  In  this 
thesis,  the  times  at  which  demands  occur  are  assumed  to  be  generated  by  independent 
Poisson  processes  with  intensities  X:  (annual  demand  rate  for  item  i).  We  assume  that 
all  shortages  are  backordered  and  will  be  satisfied  from  the  next  shipment. 

We  assume  that  a  holding  cost  of  H'  per  unit  time  is  charged  for  each  unit  of 
item  i  in  inventory.  When  an  order  is  placed,  a  setup  cost  for  an  independent 
replenishment  of  item  i  is  assumed  to  be  Cj  *  A  +  where  A  is  a  fixed  cost  and  a^  is 
a  cost  which  depends  on  the  item  ordered.  The  cost  of  a  joint  replenishment  of  both 
item  i  and  j  is  A  +  ■+•  aj  .  For  instance,  if  item  i  triggers  a  replenishment,  the  setup 

cost  of  A  +  a^  is  applied  to  item  i,  and  aj  to  the  jointly  replenished  item  j  . 

This  type  of  cost  structure  is  particularly  appropriate  when  a  group  of  items  is 
ordered  from  the  same  supplier,  or  if  a  group  of  items  uses  the  same  means  of 
transportation,  and  transportation  costs  do  not  increase  proportionally  when  the 
quantity  transported  is  increased.  Thus,  when  a  replenishment  is  made,  there  is  a 
major  fixed  cost  independent  of  the  number  of  items  involved  and  a  variable  cost  that 
depends  upon  the  number  of  items  [Ref.  1:  pp.  278-83]. 

Because  of  this  cost  structure,  it  is  clear  that  there  may  be  an  opportunity  to 
reduce  total  variable  costs  if  one  takes  advantage  of  joint  replenishment. 


B.  SCOPE 

Our  interest  in  this  problem  stems  from  the  'random  joint  order  policy'  proposed 
by  Balintly  [Ref.  2],  and  from  an  (S,cts)  policy  treated  by  Silver  [Ref.  3J.  For  the 
backorder  case  with  Poisson  demands  and  constant  lead  times,  -  and  the  specified 
customer  service  levels,  we  will  compare  the  individual  order  policy  versus  joint  order 
policy. 

However,  for  selecting  the  control  parameters  of  (S,c,s)  in  Joint  Order  Policy,  we 
have  considered  slightly  different  procedures  suggested  by  Schacck  and  Silver  [Ref.  4], 
and  Silver  [Ref.  3]. 

Two  programs  were  written  in  order  to  conduct  the  simulation  experiments.  The 
first  program  determines  the  control  parameters  of  each  item  in  the  group  for  both 
independent  and  joint  replenishment  policies.  A  second  program  evaluates  the  model. 
Detailed  user  instructions  and  the  program  itself  are  contained  in  Appendix  A  and  B, 
.while  examples  of  their  use  are  presented  in  Chapter  II. 

The  individual  order  policy  is  demonstrated  only  for  the  purpose  of  comparison. 
For  the  sake  of  simplicity,  we  compare  the  joint  order  policy  to  the  individual  order 
policy  where  the  classical  lot  size  formula  is  used  for  each  item. 


II.  GENERAL  FORMULATION  OF  THE  MODEL 


A.  ASSUMPTIONS  AND  NOTATIONS 

1.  Assumptions 

We  are  dealing  with  a  group  of  items  where  the  cost  of  replenishing  two  or 
more  items  at  the  same  time  is  less  than  the  total  cost  of  replenishing  the  same  number 
in  separate  individual  replenishments.  The  assumptions  of  the  model  are: 

1.  There  is  fixed  cost  A  associated  with  each  replenishment,  and  variable  cost  a. 
associated  with  each  item  involved  in  the  replenishment. 

2.  Demands  for  item  i  are  Poisson  distributed  with  parameter  X-. 

3.  There  are  no  quantity  discounts. 

4.  Inventory  holding  costs  are  proportional  to  the  average  dollar  value  of 
inventory. 

•  * 

5.  Service  level  is  defined  in  one  of  two  wavs:  probabilitv  of  no  shortage  per 
replenishment  cvcle,  and  fraction  of  demand  to  be  satisfied  directiv  fro“m  the 
shelf. 

2.  Basic  Notation 

We  use  the  following  notation: 

A  :  fixed  cost  per  replenishment,  (independent  of  the  number  of  items  and 
units  involved) 

aj :  fixed  cost  for  including  item  i  in  a  replenishment. 

Cj :  ordering  cost  per  order:  (  Cj  =  A  +  aj) 

ECj:  expected  relevant  annual  inventory  cost  for  item  i. 

Pj  :  desired  customer  service  level,  measured  in  terms  of  the  probability  of  no 
shortage  per  replenishment  cycle. 

P2  :  desired  customer  service  level,  measured  in  terms  of  the  fraction  of 
demand  satisfied  directly  from  stock  on  hand. 

L  :  replenishment  lead  time  in  years, 
r  :  inventory  holding  cost  rate, 
v-  :  standard  unit  price  of  item  i. 

Sj :  reorder  or  must-order  point  of  item  i. 
c-  :  can-order  point  of  item  i. 

Sj :  order-up-to  point  of  item  i. 

Qj  :  order  quantity  of  item  i. 

Xj :  Poisson  demand  rate  for  item  i,  in  pieces  per  year. 


I*  :  average  on  hand  inventory  level  of  item  i. 
n  :  number  of  items  in  the  replenishment  group. 

Nj  :  expected  number  of  replenishments  which  involve  item  i  per  year. 

NTj :  expected  number  of  replenishments  triggered  by  item  i  per  year. 

Other  notation  will  be  introduced  as  they  needed. 

B.  INVENTORY  THEORY  BACKGROUND 

Before  considering  the  coordinated  control  problem  it  might  be  helpful  to  briefly 
discuss  the  related  single  item  models  whose  solution  will  be  a  key  element  in 
determing  the  parameters  of  the  (S,  c,  s)  control  system. 

1.  Economic  Order  Quantity  Model 

The  E.O.Q.  model  is  based  on  the  assumption  that  the  entire  lot  is  added  to 
stock  at  one  time,  and  that  the  stock  will  be  withdrawn  at  a  constant  rate,  and  no 
stockouts  are  permitted. 

Total  annual  variable  cost  *  Ordering  cost  +  Holding  cost 

ECj-  XiCi/Qi  +  -Hj  Qj  /  2  (2.1) 

The  annual  order  cost  is  obtained  by  multiplying  the  number  of  orders  per 
year  (X-/Q-)  by  the  cost  of  placing  an  order  (C-). 

The  average  annual  holding  cost  for  item  i  is  the  average  inventory  (Qj/2) 
times  the  annual  unit  holding  cost  (Hj  =  vjr). 

To  obtain  the  minimum  cost  lot  size  (E.O.Q.),  take  the  first  derivative  of  total 
annual  cost  with  respect  to  Q-  and  set  it  equal  to  zero.  Solving  this  for  Qj,  we  get  the 
E.O.Q.  formula: 


Qi*  -  V  t2CiXi Hj] 


(2.2) 


Once  the  economic  order  quantity  is  known,  the  expected  number  of  orders 
placed  during  the  year,  n,  can  be  determined: 


T  V'Qi 


* 


2.  Exact  Model  for  the  Backorders  case  with  Poisson  Demands  and  Constant  Lead 
Times 


If  demands  are  random,  the  probability  cf  stockouts  exists,  and  we  must 
include  a  stockout  cost  in  the  expression  for  expected  annual  variable  costs.  This 
expression  is  given  in  [Ref.  5:  eqn.  (4-61 )]  to  be 


EC-  XiCi/Qi+  Hi  (Qj  /  2  +  1/2  +  Si  -  Mj)  +  <Pi  E(Qif  sj) 


+  (ti  +  Hi)B(Qi,  Si) 
where 

q>i :  the  cost  of  backorder  per  unit  quantity  of  backorder. 

ti :  the  cost  of  backorder  which  is  proportional  to  the  length  of  backorder  time. 

:  expected  lead  time  demand  for  item  i  in  pieces,  Mi*  XjL. 

H: :  holding  cost  of  item  i. 

*  •  •  •  • 

E(Qi,  :  expected  number  of  backorders  incurred  per  year  for  item  i. 

B(Qi,  S:) :  the  steady  state  expected  number  of  backorders  for  item  i. 

It  is  not  easy  to  derive  the  optimal  Qi  and  Si  by  use  of  the  above  exact 
expression  for  ECj.  Fortunately,  it  turns  out  in  practice  that  it  is  seldom  necessary  to 
use  the  exact  formulation,  except  for  the  case  when  it  costs  very  little  to  incur 
backorders.  The  approximation  discussed  in  the  next  section  has  been  shown  to  give 
good  results  when  backorders  are  infrequent. 

3.  Approximate  Form  with  Known  Backorder  Cost 

The  average  annual  variable  cost  equation  can  be  approximated  as  follows 
[Ref.  6:  pp.  59  -  62]: 


EC-  XiCi;Qi+Hi(Qi/2  +  si-  Mi)  +  9i(Xi/Qi)E(Mi>si) 


where 


<p*  :  the  cost  of  backorder  per  unit  quantity  of  backorder. 

E(Mj>  sp  :the  expected  number  of  backorders  at  the  end  of  cycle. 


To  minimize  ECj,  take  the  partial  derivative  of  the  ECj  with  respect  to  Qj  and 
Sj  and  set  these  equal  to  zero. 


dECj/dQ-  -  XjCj/Qj-  +  Hj/2  +  9iXiE(Mi>si)/Qi—  0 


dEq/dtj-  Hi  +  (fiAj/  Qi)dE(Mi>$i)/d$i-  0 


(2.6) 


The  first  equation(2.5)  yields  the  optimal  Qi  for  a  given  reorder  point  sj: 


Q;*  -  V  +  *iE(Mi>  Sj)  /  Hj]  (2.7) 

The  second  equation(2.6)  yields  : 

dE(Mi>si)/dsi  -  -lPr(Mi)  (2.8) 


where 

E  (Mj>  S:)  :  the  expected  number  of  backorders  at  the  end  of  cycle. 

*  Y  Pf  (M;) :  complementary  cumulative  distribution. 

This  is  the  complementary  cumulative  distribution  of  the  Poisson  random 
variable,  X,  evaluated  at  s-.  Therefore,  solving  the  second  equation  for  s-  in  terms  of 
Q-  gives  the  maximum  allowable  probability  of  a  shortage  during  a  lead  time  : 


F'(Sj)  -  (HWHvfa)  (2.9) 

To  find  the  optimal  pair  (Qj,  Sj)  that  minimize  ECj,  iterative  procedures  can  be 
applied,  such  as  described  in  [Ref.  6:  p.  61].  Hower,  we  do  not  apply  this  kind  of 
procedure  to  find  out  the  values  of  control  parameters  in  (S,  c,  s)  policy,  because  the 
order  size  in  the  (S,  c,  s)  policy  is  not  fixed. 

In  (S,  c,  s)  policy,  the  order  size  is  determined  as  an  order-up-to  point  minus  a 
reorder  point  for  the  item  which  triggers  the  replenishment  action.  Order  sizes  for 
jointly  repienishabie  items  are  order-up-to  points  minus  on-hand  levels  which  are  less 
than  or  equal  to  can-buy  points  but  more  than  reorder  points. 

4.  Some  Useful  Formulas  for  Obtaining  the  Imputed  Backorder  Cost  /  Reorder 
Point  /Safety  Quantity  in  Fixed  Order  Size  Sy  stem 

In  any  practical  situation,  it  is  very  difficult  to  determine  accurately  backorder 
costs.  They  can  include  such  factors  as  loss  of  customers'  goodwilI(i.e.,  in  the  future, 
he  may  take  his  business  elsewhere),  or  in  the  military  supply  system,  the  cost  of 
having  some  first  line  weapon  system  inoperative  because  of  lack  of  parts.  Other 
contributions  to  shortage  cost  can  be  somewhat  easier  to  measure.  However,  these 
are  usually  small  part  of  the  total  backorder  costs  [Ref.  5:  p.  18]. 
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We  shall  assume  here  that  there  is  a  backorder  cost  associated  with  each  unit 
backordered. 

A  complete  cost  minimization  analysis  of  the  optimal  values  of  the  model 
parameters  can  be  made  only  if  a  cost  can  be  assigned  to  each  shortage.  As  an 
alternative  to  specifying  the  shortage  cost,  one  can  examine  a  given  inventory  policy 
and  obtain  the  imputed  cost  of  shortage  within  it  [Ref.  7:  p.  338].  Rather  than 
explicitly  costing  backorders,  we  will  specify  a  service  level  constraint  and  determine 
the  imputed  stockout  cost. 

Whenever  an  organization  uses  a  service  level  to  establish  a  reorder  point,  it 
effectually  establishs  a  stockout  cost.  Associated  with  a  given  service  level  is  an 
imputed  or  implicit  stockout  cost.  It  is  a  simple  matter  to  determine  the  imputed 
stockout  cost  for  a  given  service  level  from  the  probability  of  a  stockout  [Ref.  8:  p. 
162]. 

a.  Probability  of  no  Shortages  per  Replenishment  Cycle 

A  service  level  based  on  the  frequency  of  service  'per  order  interval  or 
replenishment  cycle  will  indicate  the  probability  of  not  running  out  of  stock  during  the 
replenishment.  This  approach  does  not  concern  itself  with  how  large  the  shortage  is, 
but  only  with  how  often  it  can  occur  during  'he  lead  time  (replenishment  cycle).  It  is 
defined  as  the  fraction  of  the  replenishment  cycles  without  depletion  of  stock-. 

number  of  order  periods  with  a  stockout 

Service  level  fraction  per  cycle  *  l  -  - 

total  number  of  order  periods 
-l-Pr(Mi>si) 

Pr(Mi>si)  =  Pr(sp 

*  1  *  (service  level  fraction  per  cycle) 

The  term  Pr(s-)  is  the  stockout  level  fraction  per  order  cycle  or  the 
probability  of  at  least  one  stockout  while  awaiting  a  supplier's  delivery.  It  is  also  a 
measure  of  the  fraction  of  lead  time  periods  during  which  the  demand  will  exceed  the 
reorder  point.  The  magnitude  of  the  stockout  is  ignored  with  this  approach. 

The  following  formula  is  developed  to  compute  the  imputed  backorder  cost: 

Pr(Si)  -  (HjQj)  /  (bjXj)  (2.10) 


where 

Pr(sj):  probability  of  stockout  for  item  i. 
Hj :  holding  cost  of  item  i,  (H^=  v^r). 


b- :  imputed  backorder  cost  of  item  i. 

To  obtain  the  appropriate  stockout  cost,  solve  eqn.(2.10)  for  b-.  An 
example  will  illustrate  the  procedure  in  a  later  section. 

b.  Fraction  #/  Demands  Satisfied  Directly  from  the  Shelf 

The  fraction  of  units  demanded  and  immediately  filled  from  the  shelf  can  be 
defined  as  the  ratio:  (number  of  units  supplied)/(total  number  of  units  demanded) 
which  is  equal  to  {1  -  [(number  of  units  short)/(total  number  of  units  demanded)]}. 

The  above  relationships  must  be  measured  over  some  fixed  time  period 
such  as  the  order  interval. 

For  example,  to  obtain  the  stockout  level  fraction  for  units  demanded: 


Stockout  level  fraction  for  units 


E(Mi>Sj) 


quantity  demanded  during  a  cycle 


E(Mi>si)  /  Q{  -  (HjQj)  /  (bjty 


(2.11) 


where 

E(Mj>  s-):  expected  number  of  backorders  for  item  i  during  the  cycle. 

:  holding  cost  of  item  i,  (H^»  v^r). 
bj :  imputed  backorder  cost  for  item  i. 
c.  Safety  Stock 

With  backorders,  there  is  no  loss  of  sales  since  the  customer  waits  for  the 
arrival  of  the  next  order  if  stock  is  not  available.  The  expected  safety  stock  is  defined 
as: 


Safety  stock  *  £  (s  -  M)  Pr  (M)  (2.12) 

”  s  I  Pr  (M)  -  £  M  Pr  (M)  -  s  -  E(M) 

where 

E(M) :  the  expected  lead  time  demand. 

The  number  of  backorders  per  lead  time  ■  0,  if  M  <  s. 

=  M  -  s,  if  M  >  s. 


So,  the  safety  stock  is  simply  the  reorder  point  minus  the  average  lead  time 
demand.  The  reorder  point  is  determined  so  that  the  customer  service  objective  is 
satisfied. 


MATHEMATICAL  STATEMENT  OF  THE  DECISION  PROBLEM  FOR 
THE  JOINT  REPLENISHMENT  MODEL 


1.  Cost  Equation 

The  basic  cost  equation  in  the  joint  replenishment  model  is  as  follow  [Ref.  3]: 


EC—  IjV^r  +  NT- A  +  N-a^ 


(2.13) 


The  first  term  represents  the  inventory  holding  costs  and  the  other  two 
represent  the  set  up  charges  allocated  to  item  i.  The  total  expected  costs  per  year  of  the 
group  of  items  are  given  by 


EC  -  S  EC-  £  (Ijvjr  +  NTjA  +  Nft) 


(2.14) 


2.  Service  Level  Constraint 


The  service  constraint  takes  one  of  two  forms  depending  upon  the  measure  of  service 


used. 


(1)  Probability  of  no  shortage  per  replenishment  cycle 

This  is  simply  the  probability  that  demand  for  item  i  during  the  lead  time  will 
exceed  the  reorder  point. 


Pr(s)  -  Pr(Mi>si)  -  SPr(Mi) 


(2.15) 


where 


Pr(s):  probability  of  stock  out. 


Mj  :  lead  time  demand  for  item  i. 


(2)  Fraction  of  demand  satisfied  directly  from  the  shelf 

The  expected  stockout  quantity  during  the  lead  time  is  given  by 


E(M  >  Si)-  I(Mi-Si)  P_  (Mj) 


(2.16) 


I". 


where 

E(Mj>  Sj) :  expected  number  of  stock  out  for  item  i  in  pieces 
during  the  lead  time 

Mj  ;  lead  time  demand  for  item  i  . 

A  backorder  occurs  if  and  only  if  the  lead  time  demand  (which  is  Poisson  with 
parameter  Mj)  is  greater  than  the  inventory  level  of  the  item  when  the  order  is  placed. 

D.  INDIVIDUAL  ORDER  POLICY 

Under  the  the  individual  order  policy,  there  is  a  fixed  order  quantity  for  item  i 
that  is  ordered  every  time  the  reorder  point  is  reached  . 

Safety  stock  is  needed  to  protect  against  stockouts  between  the  time  the  reorder 
point  is  reached  and  the  order  received.  As  presented  in  [Ref.  3],  the  individual  order 
policy  (called  independent  control)  is  a  special  case  of  the  joint  replenishment  policy  (S. 
c,  s)  where  Cj*  s-.  In  this  case  item  i  is  ordered  only  when  the  inventory  position  hits 
its  reorder  point,  and  the  order  quantity  is  always  of  size  Q-  *  Sj  -  s-.  No  other  items 
will  be  included  in  the  replenishment,  hence  a  set  up  cost  of  A  +  a-  is  incurred  with 
each  replenishment. 

As  mentioned  in  the  introduction,  the  main  focus  of  this  thesis  is  to  compare 
joint  replenishment  with  independent  control  to  determine  the  cost  savings  that  can  be 
achieved. 

Let  us  first  consider  independent  control.  Then,  Cj  -  s-,  and  the  optimization 
problem  is: 

Minimize  ECj  =  (Sj  -  Sj  +  1)/  2  +  (Sj  -  XjL)Vjr  (2.17) 

+  (A  +  aj)Xj  /  (Sj  -  Sj) 

subject  to 

Pr(X  *  SjlXjL)  2  Pj  (2.18) 

or 


1  -  E(  X  >  SjlXjL)  Qj>  P: 


(2.101 
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where,  X  is  the  lead  time  demand  with  a  Poisson  distribution  having  mean  /.-L. 

The  left  hand  side  of  eqn.(2.18):  P^X  5  Sj|  XjL)  is  the  probability  of  no  shortage 
per  replenishment  cycle. 

In  the  left  hand  side  of  eqn.(2.19),  the  term  E(X>s-)/Q-  is  the  stockout  fraction 
of  units  demanded.  Therefore  1  -  E(X>Sj)/Qj  is  the  fraction  of  demand  satisfied 
directly  from  the  shelf. 

The  exact  solution  to  this  problem  requires  that  the  pair  (Sj,  Sj)  which  minimize 
eqn.(2.17)  subject  to  either  eqn.(2.18)  or  eqn.(2.19)  be  found  for  each  item  i.  However, 
it  have  been  observed  by  many  researchers  (see:  for  example  Silver  and  Wilson  [Ref.  9]) 
that  a  much  simplier  sequential  approach  does  almost  as  well  from  the  standpoint  of 
minimizing  costs.  This  sequential  approximation  approach  is  what  we  describe. 

1.  Applied  Solution  Procedure 

Stepl  :  Let  Qj  =  S-  -  Sj.  Find  the  value  of  Q-  which  minimizes  EC-  ignoring 
the  service  level  constraint.  Then  given  this  value  of  S-  -  s-,  the  constraint  is  used  to 
find  the  lowest  allowable  value  of  s-.  This  first  step  leads  to  the  usual  EOQ  expression 
q.  -  Sj  -  Sj  »  V  [  2(A  +  z{)\!  Hj]. 

Step2  :  Determine  the  reorder  points. 

(1) For  the  service  measure,  probability  of  no  shortage  per  replenishment  cycle, 
find  the  smallest  value  of  Sj,  which  satisfies  the  condition  of  eqn.(2.18): 

Pr(X  £  Sj  |  XjL)  2£  Pj 

(2) For  the  service  measure,  fraction  of  demands  satisfied  directly  from  the 
shelf,  find  the  smallest  value  of  sj,  which  satisfies  the  following  inequality: 

kjL  -  Sj  -  Pr(XSsj.l  |Xjl)  +  SjPr  (X  SI  Sj  |  XjL)  £  Qj(l  -  P2) 

Step3  :  Set  Sj  *  Sj  +•  Qj. 

Step4  :  Compute  the  average  annual  variable  costs. 

ECj*  (Sj  -  Sj  +  1)/  2  +  (Sj  ■  XjL)Vjr  +  (A  +  aj)Xj (Sj  -  Sj)  using  the  values 
Sj  and  Sj  estimated  from  Steps  2  and  3.  A  computer  program  for  solving  for  these 
parameters  is  provided  in  Appendix  A. 

Note  that  we  can  now  obtain  imputed  stockout  costs. 

Step5  :  Compute  the  imputed  stockout  cost 


(1) For  the  service  measure,  probability  of  no  shortage  per  replenishment  cycle: 
Use  eqn.(2.10 y.  P^s)  -  (HjQj)  /  (bjXj)  and  solve  for  bj,  bj-  HjQj  /  (1-  Pj)  Xj. 

(2) For  the  service  measure,  fraction  of  demands  satisfied  directly  from  the 
shelf,  it  is  necessary  to  solve  for  E(M  >  sj)  and  find  the  associated  P^s)  before  using 
eqn.(2.11).  It  can  be  done  alternatively  by  finding  the  complementary  cumulative 
distribution  of  the  lead  time  demand,  introduced  by  eqn.  (2.8).  Then  eqn.(2.1l)  yields 
bi*  /  [- 1  Pr  (M^)] 

2.  Numerical  Examples 

Repeated  from  Silver's  example  in  his  paper  [Ref.  3). 

Example  1.  Service  Level  =  Probability  of  no  shortage  per  replenishment 

cycle. 

Consider  the  following  example  involving  four  items: 

A  -  S  50 

a-  -  S  10  for  i  -  l,  2,  3.  4 
r  -  0.2  S  S  .  yr. 

L  -  I  month  (i.e.  1,12  year  ) 

TABLE  l 
INPUT  DATA 


Use  of  eqn.  (2.15),  (2.17),  (2.18)  gives  the  results  shown  in  Table  2. 

Example  2.  Service  Level  ■*  Fraction  of  demand  satisfied  directly  from  the  shelf. 

Consider  the  same  data  as  example  I,  but  with  a  new  the  measure  of  service 
level.  Suppose  the  specified  fraction  of  demands  to  be  satisfied  directly  is  0.99.  The 
results  shown  in  Table  3. 


TABLE  2 

INDEPENDENT  ORDER  POLICY:  SERVICE  LEVEL  95%  (Pj) 


Item 

LT  demand 

EOQ  (integer) 

Reorder 

TVC(S-) 

S/O  cost 

4 

24.167 

MB 

10.167 

m\ 

9|j99 

| 

HI 

232.024 

35.343 

88.860 

84.977 

16:162 

22.051 

13.453 

TOTAL :  441.203 


TABLE  3 

INDEPENDENT  ORDER  POLICY:  SERVICE  LEVEL  99%  (P2) 


Item 

LT  demand 

EOQ  (integer) 

Reorder 

TVC(S) 

S,  O  cost 

1 

24.167 

f 159)  ’ 

25 

220.984 

1.983 

2 

3.417 

143.178 

143) 

3 

34.383 

1.SS2 

5 

6.417 

108.840 

7 

85.740 

3.499 

4 

10.167 

178) 

10 

82.217 

1.536 

TOTAL :  423.323 


E.  JOINT  ORDER  POLICY 
1.  Introduction 

Let  us  now  consider  the  possibility  of  joint  replenishments.  As  discussed 
earlier,  we  can  reduce  ordering  setup  costs  by  combining  the  orders  of  several  items 
whenever  an  order  is  placed.  This  may  result  in  increased  holding  costs  but  the  savings 
in  ordering  costs  may  be  greater  than  the  increase  in  holding  costs  resulting  in  a  net 
decrease  in  total  costs.  Further,  it  is  likely  that  the  occurance  of  stockouts  will  be 
reduced  since  we  will  frequently  reorders  some  items  before  the  inventory  levels  hit  the 
normal  reorder  points  for  those  items. 

We  do  not  seek  to  determine  the  "optimal"  joint  replenishment  policy. 
Rather,  we  restrict  attention  to  a  simple,  intuitively  appealing,  policy  and  seek  to  find 
the  optimal  values  of  parameters  for  the  selected  policy.  The  policy  we  consider  is  the 
(S,  c,  s),  "can  buy,  must  buy"  policy  which  requires  an  order  to  be  placed  for  item  i  up 
to  the  level  any  time  the  inventory  position  for  item  i  reaches  or  falls  below  sj.  In 
addition,  however,  we  order  any  other  item  i  whose  inventory  position  is  below  the  can 
buy  point  c-v  we  always  order  up  to  the  level  S-. 


2.  General  Concepts  for  Obtaining  the  Values  of  Control  Parameters 

Suppose  a  single  item  i,  is  in  a  group  of  coordinated  items  and  assume  initially 
that  the  replenishment  lead  time  is  zero.  The  assumption  of  unit  sized  Poisson 
demands  together  with  a  zero  lead  time  implies  that  s;>  0. 

Now  item  i,  from  time  to  time,  will  be  faced  with  the  opportunity  of  a 
replenishment  at  the  reduced  setup  cost  of  only  ”aj“.  Such  an  event  is  caused  by 
another  item  triggering  a  replenishment.  Occasionally  item  i  will  trigger  an  order  when 
its  inventory  hits  the  zero  level.  Under  such  circumstances  we  would  like  to  determine 
the  (S,  c)  pair  which  minimizes  the  expected  costs  per  unit  time. 

Silver  shows  that  the  opportunities  to  replenish  at  the  reduced  cost  occur 
probabilistically  according  to  a  Poisson  process  with  a  rate  p  per  year,  where  p  is  the 
expected  number  of  orders  triggered  per  year  by  all  other  items  in  the  group. 

Let  NTj  be  the  number  of  orders  triggered  per  year  by  item  i.  Then  Pj  is  given 
by  Pj  *  H  NTj,  where  the  summation*  is  taken  over  all  i  #  j. 

For  example,  if  there  are  four  items  in  the  group  and  we  are  considering  item 
2,  we  would  have  ^  »NTj  +  NT 3  +  NT4. 

It  is  very  difficult  to  determine  the  value  of  Pj  algebraically.  We  will  consider 
some  heuristic  methods  to  approximate  this. 

Silver  [Ref.  3:  eqn.  13]  determines  that  the  relevant  cost  equation  (the 
subscripts  i  have  been  suppressed)  is  given  by: 

ECj-  (S  -  c  +  p  (1  -  pc  )  /  (1  -  p)}*1  ((S  -  cXS  +  c  +  1)  vr  /  2  (2.20) 

+  P  (c  -  p  (t  -  pc  )  /  (1  •  p)J  vr  /  (1  -  p)  +  X.  pc  A  +  X  a) 
where 

p  -  X  /  (X  +  p) 

The  algorithm  suggested  by  Silver  in  [Ref.  3]  is  a  simultaneous  optimization 
procedure  to  minimize  eqn.  (2.17)  with  respect  to  S-,  Cj,  sj.  The  procedure  is  iterative. 

Solution  procedures  to  determine  the  control  parameters  of  the  (S,  c,  s)  policy- 
are  explained  in  the  next  chapter. 


III.  DEVELOPMENT  OF  THE  MODEL 


In  this  chapter  we  consider  the  general  joint  replenishment  problem  with  positive 
leadtimes.  We  obtain  values  of  the  parameters  s-,  Cj,  and  S'  using  an  iterative 
procedure  recommended  in  [Ref.  3].  Silver  [Ref.  10]  derives  the  following  expression  for 
the  total  expected  costs  per  year  for  item  i  (we  suppress  the  item  subscript  for 
simplification  of  notation): 

Eq  -  {S  -  c  +  p  (1  .  pc)  /  (1  -  p)}’1  {(S  -  cXS  +  c  +  1)  vr  /  2  (3.1) 

+  p  [c  -  p  (1  -  pc)/(  1  -  p)]vr/(l  -  p)  +  X  pc  A  +  X  a} 

where 

p  *  X/(X  +  p) 

p  is-  the  probability  that  a  particular  occurance  is  either  a  demand  or  an 
opportunity  to  replenish  at  reduced  cost. 


A.  SOLUTION  PROCEDURE 

We  recommend  as  a  solution  procedure  a  somewhat  modified  version  of  the 
procedure  suggested  by  Silver. 

1.  Algorithm 


Step  1  :  Initialization 

Use  the  values  computed  by  the  independent  control  to  initialize.  Let  k 
where  k  *  index  of  iteration. 


1. 


cmin  “  °- 


cmax  ”  Q 


EC(k)  •  EC(determined  from  the  independent  control) 
N(k)  «  X  /  Q* 

Pj(k)  -  £  Nj,  where  all  i  *  j.  and  i,  j  =  1, ...,  n 
P(k)  -  X  /  (X  +  n) 
c(k)  -  0 
S(k)  -  Q" 
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Step  2  :  Let  k  •  k  +  1 

<*k)  -  c(k-l)  +  (Cmax  -  Cmin)  /  2.  if  EC(k)  <  EC(k-l) 
c(k)  -  c(k-l)  -  (Cmax  -  Cmin)  /  2.  if  EC(k)  >  EC(k-l) 

S(k)  -  c(k)  -  p  (1  -  pc)/(l  -  p)  +  V  {2  X  (a  +  A  Xc)/vr 
+  2  c  Xc+  !/(l  -  p)  -  p  (1  -pcXl  +  Pc+  !)/(l  -  P)2} 

EC(k)  -  {p  (1  -  pc)/(l  -  p)}*1^  [c  -  p  (1  -  pc)/(l  -  p)]vr 
/(  1  -  p)  +  k  pc  A  +  X  a} 

Step  3  :  Test  for  convergence  of  c  and  S. 

Round  the  c  and  S  values  to  the  nearest  integer. 

If  c(k)  *  c(k-l)  and  S(k)  *  S(k-I),  go  to  Step  4,  else  go  to  Step  2.  And, 
repeat  until  S  and  c  converge. 

Step  4  :  Repeat  Steps  1,  2  and  3  for  i  *  I,  2, ...,  n. 

Step  5  :  Let  Nj  -  X  p  C/{S  -  c  +  p  (1  -  p  c)/(l  -  p)}. 

p  values  are  from  the  previous  iteration  and  S  and  c  values  are  from  step  4. 

Step  6  :  Repeat  steps  1-5  until  |c(k)  -  c(k-l)|  < 6  and  |S(k)  -  S(k-l)|  < 6. 

Step  7  :  Now  determine  the  must-buy  points,  sj,  for  i  *  1 .  n  from  the  service 

constraints. 

(a)  For  Service  Measure  1:  probability  (Pj)  of  no  shortage  per  cycle,  find  the 
smallest  value  s  such  that 

(l/p)Pr(X£s  +  c|XL)-  Pr(X»s+l|X  L) 

•Ps+1Ipr<x“*0|X  LMl/p)^  2:  Pj,pc 

(b)  For  Service  Measure  2:  fraction  (P->)  of  demand  to  be  satisfied  directly 
from  shelf,  find  the  smallest  value  s  such  that 

pc  {X  L  -  s  X  L  Pr  (X  £  s  -  1  |  X  L)  +  s  Pr(X  £  s|X  L)} 

+  (1  -  P)c  +  S  z  P  *w0  {X  L  -  w0  -  X  L  Pr(w0  -  l|k  L) 

+  Pr(Xiw0)} 

*  {S-c  +p(l-pc)/(l-p)>(l-P2) 

Step  8  :  Recompute  the  values  S^  and  Cj. 

Si  "  S  +  ^ 

Cj  -  C  +  Sj 

Stop 
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2.  Flowchart 


Figure  3.1  Determination  of  coordinated  pair  (S,c). 
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IV.  SIMULATION  STUDY 


In  order  to  determine  how  costs  of  ordering,  holding,  shortage,  etc.  are  affected 
by  joint  replenishment,  a  simulation  model  was  developed.  The  purpose  of  the 
simulation  study  was  to  test  the  cost  savings  by  comparing  the  costs  for  the  two 
inventory  policies,  independent  and  joint  replenishment.  The  simulation  was  written  in 
FORTRAN77. 

A.  ASSUMPTIONS  FOR  THE  SIMULATION  STUDY 

1.  Poisson  demands  (exponential  interarrival  times). 

2.  Items  demanded  one  at  a  time. 

3.  Constant  lead  times  (independent  of  the  size  of  the  order  and  the  number  of 
items  in  a  replenishment ). 

4.  For  the  joint  replenishment  case,  place  an  order  whenever  anv  item  hits  its 
must-buv  point  s-  and  order  up  to  S;.  Also,  include  item  j  ¥>  i  if  the  inventor.- 
position  for  item)  is  below  the  can-buy  point  Cj.  Order  up  to  Sj. 

5.  For  the  independent  control,  order  up  to  level  S;  whenever  the  stock  on  hand 

lor  item  i  reaches  Sj.  1 

6.  For  a  given  item,  there  is  never  more  than  a  single  order  outstanding  in  the 
joint  replenishment  case. 

B.  MODEL  INPUTS 

1.  Number  of  items. 

2.  Demand  rate  X-  for  each  item. 

3.  Holding  cost  rate  vj  for  each  item. 

4.  Lead  time  (in  months)  L. 

5.  Type  of  service  level  constraint  (Typel.  or  Type2) 

6.  Service  level  required. 

7.  Unit  price  UP-. 

8.  Group  ordering  cost  A. 

9.  Individual  ordering  cost  aj. 

10.  Stockout  cost  rate  per  unit  stockout  (percentage  of  unit  price). 

1 1.  Time  weighted  stockout  cost  (percentage  of  unit  price). 

12.  Sj,  Cj,  Sj  for  each  item. 

13.  Model  options  (Independent  control  or  Joint  Replenishment,  etc.) 

14.  Limits  of  simulation  period,  NSI.M. 


C.  MODEL  OUTPUTS 

1.  Achieved  service  level  for  each  item. 

2.  Number  of  orders  (independent  and  joint). 

3.  Total  ordering  cost  for  each  item. 

4.  Total  holding  cost  for  each  item. 

5.  Total  unit  years  of  backorders. 

6.  Total  number  of  backorders. 

7.  Total  unit  years  of  stock  held. 

8.  Total  stockout  cost  (fixed  and  time  weighted)  for  each  item. 

9.  Total  annual  costs  for  each  item. 

10.  Total  annual  variable  costs  for  all  items. 

11.  Standard  deviation  of  total  annual  variable  costs. 

12.  A  complete  audit  trail  of  all  orders  placed  (optional). 

D.  STRUCTURE  OF  THE  PROGRAM 
1.  Algorithm 

Step  1  :  Read  input  data,  process  options  and  initialize  variables. 

Step  2  :  Determine  the  earliest  event  time  and  type  of  event. 

Let  TDAj  be  the  time  of  the  next  demand  for  item  i  and  TO  A-  be  the  next 
time  of  the  for  item  of  an  order  for  item  i. 

Set  the  master  clock  time,  CT,  to  the  minimum  of  {TDAj,  TOA-). 

If  the  event  type  is  a  demand  for  item  i  then  generate  another  demand  time 
for  that  item. 

If  the  event  type  is  an  order  arrival  then  increament  the  on  hand  inventory  for 
all  items  included  in  the  order. 

Step  3  :  Determine  if  an  order  is  required. 

If  the  onhand  level  is  less  than  or  equal  to  the  reorder  point  for  that  item  and 
there  are  no  outstanding  orders  for  that  item,  place  an  order.  If  the  ordering  policy  is 
independent,  go  to  step  5. 

Step  4  :  Test  for  joint  replenishments. 

If  the  onhand  level  is  less  than  or  equal  to  the  can  order  point  for  another 
item  j  in  the  replenishment  set,  place  joint  replenishment  for  the  items  to  the  level  S-. 
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Step  5  :  Determine  if  the  simulation  is  complete. 

IfCT  £  NSIM  stop,  else  go  to  step  2. 

The  program  collects  detailed  statistics  at  each  event  time  for  unit  years  of 
stock  held,  number  of  stockouts,  unit  years  of  stockouts,  number  of  orders,  number  of 
units  ordered,  and  number  of  demands  for  each  item. 

At  the  conclusion  of  the  simulation  additional  summary  statistics  are 
computed  for  achieved  service  levels  and  costs. 


2.  Flowchart 


Figure  4. 1  Simulation  Flowchart. 


V.  EVALUATION  OF  THE  MODEL 


A.  PERFORMANCE  ANALYSIS 
1.  Basic  Examples  Tested 

Table  4  shows  the  data  used  for  the  evaluation  of  the  joint  replenishment 
model.  Several  combinations  (96  cases)  of  service  level,  service  measure,  and  lead  times 
were  considered.  The  data  include  the  test  set  considered  in  Silver's  paper  [Ref.  3]. 

TABLE  4 

ITEM  CHARACTERISTICS 


group 

number 

A 

(S) 

a 

(S) 

number 
|  of  items 

i 

,  ITEM  DATA 

A.;  (piece/yr.)  v;  (S/piece) 

1 

50 

10 

6.90 

1.20 

3.90 

2.30 

2 

125 

80 

10 

I 

40.630 

0.45 

2 

4.090 

36.41 

3 

34.680 

42.17 

4 

4.240 

4.46 

5 

4.240 

4.46 

6 

4.090 

36.41 

7 

4.240 

4.46 

8 

28.780 

3.92 

9 

4.090 

36.41 

10 

4.090 

36.41 

3 

50 

5 

4 

l 

290.000 

6.90 

-> 

4v 

41.000 

1.20 

3 

77.000 

3.90 

4 

122.000 

2.30 

4 

50 

5 

8 

1 

290.000 

6.90 

2 

41.000 

1.20 

3 

77.000 

3.90 

4 

122.000 

2.30 

5 

50.000 

1.20 

6 

154.000 

3.90 

7 

87.000 

2.30 

8 

25.000 

1.20 

2.  Model  Results 

a .  Convergence  of  Parameter  Values  and  Costs 

Table  5  shows  the  results  of  each  iteration  from  Silver's  model  for  joint 
replenishment.  Observe  that  five  iterations  were  required  for  convergence.  This 


example  considers  service  measure  type  2  with  a  specified  level  0.99.  The  lead  time  is 
one  month. 

TABLE  5 

PARAMETERS  AND  COSTS 


b.  Cost  Savings:  Joint  Replenishment  VS.  Independent  Control 

Table  6  shows  the  cost  savings  achieved  by  the  joint  replenishment  policy 
for  both  the  modified  Silver's  method  and  our  heuristic  method. 

The  comparisions  are  provided  for  two  cases: 

1.  Service  measure  1,  Pj  ■  0.95 

2.  Service  measure  2,  P2  *  0.99 

Observe  that  both  the  heuristic  method  and  Silver's  method  achieve  cost 
savings  greater  than  13  %  by  using  joint  replenishment,  and  the  heuristic  method 
results  in  greater  cost  savings  than  does  Silver's  method. 

3.  Simulation  Results 

Figure  5.1  shows  the  time  history  of  stock  on  hand  for  items  1  to  4  generated 
by  one  simulation.  The  upper  graph  represents  individual  control  and  the  lower  graph 
represents  joint  replenishment. 


TABLE  6 

COMPARISON  OF  TWO  METHODS 


1.  modified  Silver's  method 


item 

1 

2 


3 


52 


83 


total  cost 


=  0.95,  L  =  1  month 


EC( IP) 


232. 024 
35.  243 


83. 860 


67.1722  |  84.977 


jst  !  383.0190  I  441.203 


s 

EC( JP ) 

184 

221.  6339 

100 

23. 2671 

97 

70. 9459 

156 

67. 1722 

2.  heuristic  method 


item 


=  0.95,  L  =  1  month 


Save(  %) 


4.  48 


34.  17 


20.  16 


20.  95 


I  13. 19 


P 

s 

C 

S 

EC( JP) 

EC(IP) 

Save(  %) 

0.  9942 

31 

110 

177 

211.  2798 

232.024 

8.  94 

0.  9272 

3 

34 

93 

21.  6035 

35.343 

38.  87 

0. 9650 

9 

36 

94 

68. 7113  i 

88.  860 

0. 9774 

14 

59 

152 

65.3723 

84.  977 

23.07 

total  cost  I  351.8467 


441.203  I  16.83 


3.  modified  Silver's  method 


item 


1 

— 


3 


=0.99,  L  =  1  month 


4.  heuristic  method 


item 


1  0. 9942 

2  0. 9272 
0. 9650 
0.  9774 


S 

EC( JP ) 

EC( IP) 

177 

211.  9739 

220. 984 

95 

22. 0671 

34.  383 

93 

67. 8259 

85. 740 

150 

64.  4122 

82. 217 

48 


77 


total  cost 


0.  99, 


366. 2788 


1  month 


423. 323 


Save(%) 


4.  08 


35.  82 


20.  99 


21.  66 


13.  48 


s 

C 

25 

104 

-1 

30 

5 

32 

8 

53 

S 

EC( 

JP) 

EC(  IP) 

Save( %) 

171 

202. 

5998 

220.  984 

3. 

14 

89 

20. 

6435 

34.  383 

39. 

96 

90 

65. 

5912 

85. 740 

23. 

50 

146 

62. 

6123 

82. 217 

23. 

85 

total  cost  |  351. 8467 


423. 323 


16.  88 


i 


The  simulation  results  for  several  combinations  of  leadtime,  serv  ice  measure 
and  service  level  are  represented  in  Tables  7,8,9. 

When  lead  time  is  less  than  1  month,  the  heuristic  method  yields  more  savings 
in  total  cost  and  also  satisfies  the  required  service  level.  But,  as  the  lead  time  becomes 
larger,  it  sometimes  fails  to  meet  the  required  service  level,  although  it  does  continue  to 
generate  greater  cost  savings  than  Silver's  method  with  the  cost  parameters  examined 
in  this  thesis.  Both  methods  generate  substantial  savings  over  independent  control. 
In  short  lead  times,  both  methods  tend  to  overprotect  against  stockouts.  They 
generally  produce  service  levels  in  excess  of  what  is  required. 

One  explanation  for  the  better  performance  of  the  joint  replenishment  models 
for  the  short  lead  time  examples  is  the  fact  that  a  key  step  in  the  joint  replenishment 
algorithm  for  determining  the  S,  c,  s  parameters  assumes  that  lead  time  is  zero. 
Adjustments  are  then  made  to  those  values  to  account  for  the  nonzero  lead  times  but, 
apparently  those  adjustments  are  not  adequate.  Additional  work  is  required  to 
examine  this  issue  and  to  explore  other  methods  of  adjustment. 

We  should  keep  in  mind  the  fact  that  the  model  satisfies  all  requirements  in 
the  experiment  does  not  mean  that  it  is  a  correct  model  -  merely  that  it  is  a  plausible 
one  which  has  not  been  found  inadequate  by  the  data  or  experiment. 


TABLE  7 

COST  SAVINGS  BY  JOINT  REPLENISHMENT:  WHEN  L  -  t  MONTH 


1.  PARAMETERS  FROM  MODIFIED  SILVER'S  METHOD 


Group 

number 

req. 

svS 

(%) 

AVG.  Total  Variable  COsts($) 

INDEPENDENT  JOINT 

Cost  Savings! %) 
by  Joint 
Replenishment 

AVG.  SVL 
achieved 

hv 

MODEL 

SIM 

MODEL 

SIM 

MODEL 

SIM 

SIM. 

1 

Pl=95 

441.  20 

443. 04 

383. 02 

361.  48 

13.  19 

18.  41 

99.  49 

Pl=99 

448. 06 

445. 77 

392. 74 

368.  70 

12.  35 

17.  29 

99.  74 

P2=95 

403.20 

414.  71 

343. 31 

341.  12 

14.  86 

17.  75 

97.  91 

P2=99 

|  423.32 

431.  00 

366.  28 

348.  38 

13.  48 

19.  17 

99.  71 

2 

P  1=95 

1111. 43 

1073. 85 

998.32 

953. 78 

10.  18 

11.  22 

99.  27 

Pl=99 

1153. 33 

1107.  24 

1037. 70 

956.  40 

10.  03 

13.  62 

100. 00 

P2=95 

1032. 91 

1029. 70 

878. 37 

888. 82 

14.  96 

13.  68 

97.  70 

1  P2=99 

1088. 91 

1066. 90 

971.  82 

926.  96 

10.75 

13.  12 

99.  64 

3 

P1=9S 

423.30 

422. 05 

356.  89 

345. 93 

15.  69 

18.  04 

99.  76 

Pl=99 

430. 16 

42S. 74 

365. 83 

354.  15 

14.  95 

16.  82 

100. 00 

P2=9S 

386. 08 

396.  72 

315.  71 

315.  35 

18.  23 

20.  51 

98.  44 

P2=99 

406.  80 

409. 84 

340.  01 

333. 05 

16.  42 

18.  74 

99.  76 

4 

P  1=95 

676.  10 

677. 55 

513. 14 

490. 92 

24.  10 

27.  54 

99.  21 

P  1=99 

686. 94 

683. 89 

528.  48 

508.  80 

23.07 

25.  60 

99.  80 

P2=95 

620. 56 

627. 49 

458.  34 

433. 05 

26.  14 

30.  99 

99.  01 

P2=99 

650. 70 

650. 80 

490.  60 

459. 33 

24.  60 

29.  42 

99.  83 

2.  PARAMETERS  FROM  l 

JEURISTIC 

I  METHOD 

1 

Pl=95 

441. 20 

443. 04 

366.  97 

365.  30 

16.83 

17.  55 

96.  51 

Pl=99 

448.  06 

445. 77 

376. 93 

368.  73 

15.98 

17.  28 

99.  04 

P2=95 

403. 20 

414.  71 

329. 23 

329.  36 

13.35 

20.  58 

97.  34 

P2=99 

423. 32 

431.  02 

351.  85 

347.  15 

16.88 

19.  46 

99.  46 

2 

Pl=95 

1111. 43 ! 

1073. 85 

949. IS 

960. 01 

14.  60 

10.  60 

9".  43 

Pl=99 

1153. 33 

1107. 24 

999. 63 

979. 11 

13.32 

11.  57 

99.  44 

P2=95 

1032. 91 

1029. 70 

849. 28 

845. 07 

17.78 

17.  93 

97.  54 

P2=99 

1088. 91 

1066. 90 

903. 69 

893. 75 

17.01 

16.  23 

99.  03 

3 

Pl=95 

423. 30 

422. 05 

332. 99 

342. 74 

21.33 

13.  79 

99.  60 

?1=99 

430. 16 

425. 74 

343.  17 

350. 13 

20.  22 

17.  76 

99.  77 

P2=95 

386. 08 

396.  72 

296. 95 

314.  70 

23. 09 

20.  67 

98.  01 

P2=99 

406. 30 

409. 34 

319. 11 

329.  97 

21.  56 

19.  76 

99.  64 

4 

Pl=95 

676. 10 

677. 55 

447. 50 

455. 68 

33.  81 

32.  75 

97.  19 

P  1=99 

636.  94 

683. 89 

463. 88 

471.  45 

32.  47 

31.  06 

99.  51 

P2=95 

620. 56 

627. 49 

392. 74 

394.  17 

36.  71 

37.  13 

97.  54 

P2=99 

650. 70 

650. 80 | 

|  426.24 

420.  22 

34.  50 

35.  43 

99.  29 

TABLE  8 

COST  SAVINGS  BY  JOINT  REPLENISHMENT:  WHEN  L  -  3  MONTH 


1.  PARAMETERS  FROM  MODIFIED  SILVER'S  METHOD 

Group 

number 

Wt 

(X) 

AVG.  Total  Variable  COsts($) 

Cost  Savings(%) 
by  Joint 
Replenishment 

AVG. SVL 
acLllved 
S& 

INDEPENDENT 

JOINT 

MODEL 

SIM 

MODEL 

SIM 

MODEL 

SIM 

1 

Pl=95 

Pl=99 

453. 70 
466. 64 

451.  55 
464.  38 

394.  04 
408. 46 

368.  34 
383. 34 

13.  15 
12.  47 

18.  43 

17.  45 

98.  97 

99.  74 

P2=95 

P2=99 

404.  20 
431.  48 

407. 47 
433. 39 

340.  27 
373. 06 

331.  60 
350. 99 

15.  82 
13.  54 

18.  62 

19.  01 

96.  91 

99.  59 

2 

Pl=95 

?1=99 

1173. 94 
1224.  36 

1135. 01 
1177. 53 

1019. 79 
1078. 73 

938. 40 
994.  54 

13.  13 
11.  89 

17.  32 

15.  54 

98.  33 

99.  33 

P2=95 

P2=99 

1047. 78 
1109. 34 

1054. 40 
1097. 29 

893. 13 
992. 50 

855. 98 
917.  02 

14.  76 
10.  53 

18.  82 

16.  42 

98.  34 

99.  78 

3 

Pl=95 

Pl=99 

435. 80 
448. 74 

439. 10 
446. 89 

365.  75 
382. 33 

358. 36 
371.  04 

16.  07 
14.  80 

18.  39 

16.  97 

96.  73 

99.  51 

P2=95 

P2=99 

386.  76 
413. 82 

390. 60 
410. 99 

314.  71 
347. 03 

309.  32 
337. 40 

18.  63 
16.  14 

20.  31 

17.  91 

97.  38 

99.  57 

4 

P1=9S 

Pl=99 

694.  66 
713. 84 

690. 92 
711.  89 

525.  31 
549. 89 

491.  31 
516.  17 

24.  38 
22.  97 

23.  32 

27.  49 

99.  32 

99.  80 

P2=95 

?2=99 

619. 70 
660. 82 

623. 99 
660. 38 

446. 77 
498. 85 

418. 44 
467. 03 

27.  91 
24.  51 

32.  94 

29.  28 

98.  70 

99.  84 

2.  PARAMETERS  FROM  HEURISTIC  METHOD 

1 

P1=9S 

Pl=99 

453. 70 
466. 64 

451.  55 
464.  38 

377. 98 
392. 64 

363. 73 
378.  46 

16.  69 
15.  86 

19.  45 

18.  50 

97.  35 

99.  76 

P2=95 

P2=99 

404.  20 
431.  48 

407. 47 
433. 39 

326. 78 
357.  32 

318.  15 
346.  98 

19.  15 
17.  19 

21.  92 

19.  94 

96.  75 

99.  44 

2 

Pl=95 

?I=99 

1173. 94 
1224.  36 

1135. 01 
1177. 53 

978. 83 
1041. 41 

966. 47 
1015. 42 

16.  62 
14.  94 

14.  85 

13.  77 

*93. 78 

99.  17 

P2=95 

?2=99 

1047. 78 
1109. 34 

1054. 40 
1097.  29 

854.  13 
950.  04 

360. 07 
923. 59 

13.  48 

14.  36 

18.  43 

15.  33 

97.  94 

99.  44 

3 

PI =9  5 
?1=99 

435.  30 
448. 74 

439.  10 
446. 39 

340.  91 
353. 43 

344.  09 
360. 38 

21.  77 
20.  13 

21.  64 

19.  25 

99.  20 

99.  77 

P2=95 

P2=99 

386. 76 
413- 82 

390. 60 
410. 99 

291.  57 
324.  41 

301.  49 
327. 00 

24.  61 
21.  61 

^  2  3  ? 

26: 44 

97.  70 

99.  74 

4 

Pl=95 

Pl=99 

694.  66 
713. 84 

690. 92 
711.  89 

456. 20 
483. 56 

447. 29 
474.  04 

34.  33 
32.  26 

35.  26 

33.  41 

95.  34 

98.  97 

P2=95 

P2=99 

619. 70 
660. 32 

623. 99 
660. 38 

381.  62 
432. 52 

400. 09 
425. 76 

38.  42 
34.  55 

35.  88 

35.  53 

95.  55 

99.  10 

35 


TABLE  9 

COST  SAVINGS  BY  JOINT  REPLENISHMENT:  WHEN  L  -  6  MONTH 


1.  PARAMETERS  FROM  MODIFIED  SILVER'S  METHOD 

Group 

number 

reg. 

SVC 

<%) 

AVG.  Total  Variable  COsts($) 

Cost  Savings(%) 
by  Joint 
Repleni shment 

a \rci  SVL 
achieved 
S& 

INDEPENDENT 

JOINT 

MODEL 

SIM 

MODEL 

SIM 

MODEL 

SIM 

* 

Pl=95 

?1=99 

465. 47 
485. 50 

461.92 
482. 24 

404. 32 
425. 84 

374.  87 
396.  24 

13.  14 
12.  29 

18.  85 

17.  83 

97.  53 

99.  24 

P2=95 

P2=99 

406.  25 
441. 41 

407. 04 
440.  34 

343.  22 
381.  50 

320.  37 
354.  81 

15.  51 
13.  57 

21.  29 

19.  42 

96.  23 

99.  34 

2 

Pl=95 

Pl=99 

1217. 73 
1285. 89 

1174.  03 
1240. 05 

1068. 46 
1131.  33 

983. 25 
1045. 79 

12.  26 
11.  98 

16.  25 

15.  67 

97.  90 

99.  33 

P2=9S 

P2=99 

1046. 00 
1161. 47 

1056. 59 
1132. 14 

897. 43 
1014.  36 

351.  39 
933. 91 

14.  20 
12.  67 

19.  42 

17.  51 

98.  07 

99.  77 

3 

P1=9S 

Pl=99 

447. 56 
467.  60 

442. 91 
464.  06 

375. 58 
397. 10 

359. 91 
380.  81 

16.  08 
15.  08 

18.  74 

17.  94 

94.  80 

99.  14 

P2=95 

P2=99 

390. 96 
423. 50 

392. 25 
418. 99 

317. 44 
355. 48 

310.  39 
342. 46 

18.  31 
16.  06 

20.  87 

18.  27 

97.  19 

99.  29 

4 

P  1=95 
Pl=99 

712. 21 
740. 67 

711. 50 
737. 77 

536.  58 
571.  58 

505. 25 
529. 95 

24.  66 
22.  83 

28.  99 

28.  17 

96.  47 

99.  04 

P2=95 

P2=99 

624.  71 
674.  81 

629. 18 
674.  72 

447. 04 
507. 72 

411.  24 
465. 69 

28.  44 
24.  76 

34.  64 

30.  98 

97.  19 

99.  44 

2.  PARAMETERS  FROM  HEURISTIC  METHOD 

1 

P1=9S 
P  1=99 

465. 47 
485. 50 

461.  92 
482. 24 

387. 43 
409. 79 

379. 63 
401.  59 

16.  77 
15.  60 

17.  82 

16.  72 

*92. 64 

98.  74 

P2=95 

P2=99 

406.  25 
441. 41 

407.  04 
440.  34 

328.  41 
366.  23 

328.  97 
363. 02 

19.  16 
17.  03 

19.  18 

17.  56 

95.  33 

99.  08 

2 

Pl=95 

Pl=99 

1217. 73 
1285. 39 

1174.  03 
1240. 05 

1030. 96 

1090. 96 

991.  45 
1049. 95 

15.  34 
15.  16 

15.  55 

15.  33 

98.  92 
ICO. 00 

P2=95 

P2=99 

1046. 00 
1161.  47 

1056. 59 
1132. 14 

848. 43 
962. 32 

978. 96 
940. 00 

18.  89 
17.  15 

16.  81 

16.  97 

96.  53 

59.  34 

3 

Pl=95 
P  1=99 

447. 56 
467. 60 

442. 91 
464.  06 

351.  19 
372- 95 

341.  63 
362. 02 

21.  53 
20.  24 

O'? 

4*  .  O  ' 

21.  99 

*91.  90 
*96. 14 

P2=95 

P2=99 

390. 96 
423. 50 

392. 25 
418. 99 

292. 41 
331.  47 

298.  23 
330.  47 

25.  21 
21.  73 

23.  97 

21.  13 

55.  57 

98.  48 

4 

Pl=95 
P  1=99 

712.  21 
740. 67 

711. 50 
737. 77 

466. 69 
503. 17 

532. 09 
559.  74 

34.  47 
32.  07 

25.  21 

24.  13 

*86. 62 
*94.  34 

P2=95 

P2=99 

624.  71 
674.  81 

629. 13 
674.  72 

381.  93 
442. 41 

404.  55 
452. 60 

mm 

1  ■ 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

For  each  of  several  test  combinations,  we  determined  the  independent  and 
coordinated  control  parameter  values  and  the  associated  expected  costs.  We  simulated 
several  years  of  operations  with  the  derived  joint  replenishment  ordering  parameters  to 
evaluate  the  effectiveness  of  the  models  as  compared  to  independent  control  models. 
We  also  use  the  simulation  to  compare  a  simple  heuristic  joint  replenishment  model. 
With  the  simulation  we  were  able  to  accumulate  costs  (ordering,  holding  and  stockout) 
and  effectiveness  as  measured  by  two  service  levels. 

In  summary,  the  following  statements  can  be  made: 

1.  Coordinated  control  results  in  total  costs  significantly  lower  than  under 

•  independent  control.  The  average  cost  savings  are  about  f8.7  %  . 

2.  Coordinated  control  provides  better  service  than  independent  control.  In  a  few 
cases,  the  actual  service  for  an  individual  item  is  slightly  lower  than  required: 
when  this  occurs,  the  item  involved  is  the  one  that  triggers  the  majority  of  the 
replenishments  in  the  group. 

3.  Cost  savings  of  coordinated  control  increases  as  the  ratio  a/ A  decreas  es. 

4.  Cost  savings  improve  as  the  number  of  items  in  replenishment  set  increases:  the 
more  items  there  are  in  the  group,  the  more  beneficial  coordinated 
replenishment  becomes. 

5.  For  a  given  service  measure  and  lead  time,  the  cost  savings  increase  as  the 
desired  service  level  decreases. 

B.  RECOMMENDATIONS 

The  procedure  to  determine  the  parameters  for  coordinated  control  requires  only 
the  data  for  independent  control.  We  observe  that  the  joint  replenishment  models 
perform  best  for  shorter  lead  times.  This  is  probably  the  result  of  the  fact  that  the 
joint  replenishment  algorithms  initially  assume  zero  lead  times  and  then  adjust  the 
derived  parameters  heuristically  to  account  for  nonzero  lead  times.  This  adjustment 
procedure  may  be  inadequate.  Additional  work  is  needed  to  consider  other  adjustment 
methods. 
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APPENDIX  A 

*PARA':PROGRAM  TO  SELECT  PARAMETERS 

**  VARIABLE  DEFINITION  ** 


DM 

MR 

S/Sl 

C/Cl 

IPC 

IPS 

?CI 

EC/EC1 

TLM 

TLT 

TL 

DLT 

NIT 

TN 

R  / 

F 

UP 
OCI 
OCG 
SVL 
SCU 
PSO 
STON 


RHO 


ANNUAL  DEMAND 
REORDER  POINT 
ORDER  -  UP  POINT 
CAN  ORDER  POINT 
OPTIMAL  PAIR  OF  C 
OPTIMAL  PAIR  OF  S 
ECNOMIC  ORDER  QUANTITY 

AVERAGE  ANNUAL  TOTAL  VARIABLE  COST ( INDEPENDENT  ODER) 
AVERAGE  ANNUAL  TOTAL  VARIABLE  COST (JOINT  ORDER) 

LEAD  TIME (IN  MONTH) 

LEAD  TIME (IN  YEARS) 

AVERAGE  LEAD  TIME (IN  YEAR) 

LEAD  TIME  DEMAND 

NO.  OF  TOTAL  ITEM  IN  THE  GROUP 

NO.  OF  TRIGGERING  AN  ORDER  BY  ITEM  I  IN  THE  GROUP 
OPPORTUNITY  TO  REPLENISH  AT  REDUCED  COST 
HOLDING  COST  RATE  TO  PROCUREMENT  COST 
UNIT  PRICE 

INDIVIDUAL  ORDERING  COST 
GROUP  ORDERING  COST 
SERVICE  LEVEL 

IMPUTED  UNIT  STOCKOUT  COST  GIVEN  SVL. 

MAX.  ALLOWABLE  ROBABILITY  OF  STOCK  OUT 
EXPECTED  NUMBER  OF  STOCKOUT  DURING  A  ORDER  CYCLE 


DIMENSION  ITM(IQO) , DM (100) , UP (100) ,H(100) ,0CI(100) ,ECI(100) ,Q(L00) 
&,C(100) ,TN(10Q) . S(100) , IS(100> ,R(100) ,U(100) .EC (100) , MR (100) , 
&IQ(100),DLT(100),SCU(100),IC(l00),RHo(50,100},ECl(50,100), 
&Cl(50,100) ,31(50,100) ,IC1( 50,100) ,151(50,100; , IPC (100) , IPS (100) , 
&MR1(100) , TL(lOO) 

*  DATA  INITIALIZATION 
STN=0 . 

TVC=0. 


*  READ  PROCESS  OPTIONS 

CALL  OPTION ( ID , LT , KEY , LTT , SVL , MSRCH ) 

READ (ID , 110 )NIT ,OCG, F 
TLT=LT/12. 

DO  11  1=1 , NIT 

READ(ID,111)ITM(I) ,DM(I),UP(I) , OCI(I) 
H(I)=UP(I)*F 

if(ltt.eq.i)tl1i)=tlt 

IF (LTT. EQ. 2) THEN 
CALL  LEAD ( TLT, TLY) 

TL(I)=TLY 
END  IF 

11  CONTINUE 

*  PRINT  SELECTED  OPTIONS 

WRITE (7, 100) KEY, SVL, LTT, ID, MSRCH 

*  CALCULATE  E.O.Q 

DO  1  I  =1 ,NIT 
DLT(I)=DM(I)*TL(I) 


38 


non  nnnn  n  nnnnnn 


OC*OCG+OCI ( I ) 


COMPUTE  REORDER  POINT  t IND .  ORDER  POLICY 


1  j  PROS. OF  NO  SHORTAGE  PER  REPLENISHMENT  CYCLE 

2  i  FRACTION  OF  DEMAND  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF 


CALL  POISON (KEY, DLT( I ) ,Q(I) ,SVL,MX,M,MAX,CCDF1) 

MR1 (I ) =M 
I F ( KEY . EQ . 1 ) THEN 

SCU ( I ) =H ( I ) *Q ( I ) / ( ( 1 . -  SVL)*DM(I) ) 

END  IF 

IF(KEY.EQ.2)THEN 

SCU(I)=H(I )*Q(I)/ (CCDF1*DM(I) ) 

END  IF 

ECI(I )“H(I)*( (0(I)+1 . )/2 .+(MR1 (I)-DLT(I) ) )+OC*(DM(I)/Q(I) )  / 

WRITE (7 , 501 ) ITM(I ) ,DLT(I) ,TL(I),Q(I),IQ(I) ,MR1 (I ) ,ECI (I ) ,SCU(I ) 
TVC=TVC+ECI ( I ) 

TN(I)=DM(I)/Q(I) 

STN=STN+TN ( I ; 

CONTINUE 


WRITE(7,502)TVC 


DETERMINE  THE  PARAMETERS  FOR  JOINT  ORDER  POLICY 
(REORDER  /  CAN  ORDER  /  ORDER- UP  POINT) 


DO  2  1=1, NIT 
U(I)=STN-TN(I) 

DO  3  1=1, NIT 
R(I)=DM(I)/(DM(I)+U(I)) 


*  SEARCH  LOOP  FOR  DETERMINIG  A  PAIR  (S,C)  :  ITEM  BY  ITEM 


II*l 
SEC=0 . 

SEC1=Q . 

DO  4  1=1, NIT 

RHO(I,II)=R(I) 
K=1 

CMIN=0 . 

CMAX=Q ( I ) 
C(K)=0. 
S(K)=Q(I) 
EC(K)=ECI(I) 
TEST=10.**10 
ICASE=1 
LOOP=0 
K=K+1 

IF( ICASE . EQ , 1 )TH1 


IF(ICASE.EQ.1)THEN 
C(K)=C(K-l)+(CMAX-CMIN)/2 
END  IF 

IF(ICASE.EQ.2)THEN 
C(K)=C(K-1 ) - (CMAX-CMIN)/2 
END  IF 

S(K)  =  C(K)-R(I )*( 1-R(I }**C(K) )/(l-R(I) )+SQRT(2*DM( I )*(OCI (I ) 

&  +OCG*R(I)**C(K))/H(I)+2*C(K)*R(I)**(C(K)+l)/  l-R(I)) 

&  -R<I)*U-R(I)**C(K))*(1+R(I)**(C(K)+1)>/(1-R<I))**2) 

J?  I K  ic«f  *  iMUtjtoi  5»}((i5:5(I)i  L»(I 
&/(1.-R(I)))*H(I)/(1.-R(I))+DM(I)*R(I)**C(K)*0CG+DM(I)*0CI( 
IS(K)=lNT(S(K)+.5) 

IC  K)=INT(C(K)+.5) 

IF(C(K) .EQ.C(K-l) .AND.S(K) .EQ.S(K-l) )GO  TO  6 
IF(EC(K) .GT. (0.9995*EC(K-1)) .AND.EC(K) .LE.EC(K-l))  GO  TO  6 
IF(EC(K; .LT.EC(K-l) )  THEN 

IF(C(K) .GT.C(K-1))THEN 


(S(K)-C(K))*(S«)+C(K)+lJ*H(I)/2+R(i)*(C (K)-R(I * ( 1 . -R(I )**C(K) 
l.-R(I)))*H(I)/(l.-R(f))+DM(I)*R(I)**C(K)*OCG+DM(I)*OCI(I)) 


ELSE 


CMXN-C(K) 
X CASE-1 
GO  TO  5 
ELSE  ,  v 
CMAX*C(R) 
I CASE-2 
GO  TO  5 
END  IF 


IF(C(K).LE.C<K-1))THEN 

CMIN-C(K) 

I CASE-1 
GO  TO  5 
ELSE 

CMAX-C(K) 

I CASE-2 
GO  TO  5 
END  IF 

END  IF 

IF ( LOOP. EQ.O) THEN 
CMIN-0 
CMAX-C (2) 

I CASE-2 
LOOP-1 
GO  TO  5 
END  IF 

*  FIND  OUT  THE  PAIR  (SfC),  MINIMIZES  EC 
DO  7  J=1,K 

IF (EC( J) .LE. TEST) THEN 
TEST-EC(J) 

EC1  (I , II )-EC( J) 

cki,ii)=c?j) 

SI (I, II)=S( J) 

END  IF 
CONTINUE 


ICl(I,II)=INT(Cl(I,II)+.5) 

IS1(I , II j=INT(Sl(I , II)+. 5) 

SEC1-SEC1+EC1 (I , II ) 

CONTINUE 

SHAVE- (TVC-SEC1 ) /TVC 

WRITE(7, 507)11 
DO  20  1=1, NIT 

HRITE(7 , 508)1 ,RHO(I , II ) , IC1(I,II) ,IS1(I ,11) ,EC1(I,II) 

CONTINUE 

WRITE ( 7 , 509 ) SEC1 , SHAVE 

IF(MSRCH.EQ.2)  GO  TO  42 

*  COMPUTATIONS  NEEDED  ONLY  FOR  S.  METHOD  1 

IF(II.GE.2)THEN 
DO  21  1=1, NIT 

IF (Cl (I, II). EQ. Cl (1,11-1) .AND .SI (I , II).EQ.S1(I,II~1))  GO  TO  42 
END  IF 

IF(II.GE.2)THEN 
DO  22  1=1, NIT 

IF(II .GT .NIT.OR.EC1 (I , II ) .GT. (0 . 9995*EC1 (1,11-1)) 

&  .AND .EC1 ( I , II) . LE .EC1 ( I ,11-1))  GO  TO  42 

END  IF 

IF  (II.GE.7)  GO  TO  42 
DO  23  1=1, NIT 

TN(I)  =  DM( I )  *  R(I)**C1(I.II)  /  (  S1(I. II)  -  Cl(KII)  +  R(I) 
&  *(1  -  R(I)**Cl(I,Ilj)  /  (1  -R(I)j) 

CONTINUE 

SUM=0. 


non 


dc  24  i»ijrrr 

SUK-SUIHTH  X) 

D0  25  X-lJFt 
Om-SW-TIT  X) 

do  2«  i-i.ict 

S|X)-M^X)/<DN(X)4U(X)) 

S1C1-0. 

00  TO  41 

WRITI(7,503) 

DO  27  1*1 ,NIT 
C(I  -C1?I,II) 

S  I  -SI  I, II 


Cl  -Cl  I, II) 
S(I)-S1(I,II) 
ipc(i}-:nt?c(i)+.5) 
lPS(I)«INT(S(I)+.5) 
WRITE(7,504)I,IPC(I), 
WRITS (7 ,200) 


IPS(I) 


*  COMPUTE  REORDER  POINTs JOINT  ORDER  POLICY 


DO  28  1*1, NIT 
R(I)«RHO(I.II) 
IF(KEY.EQ.1)THEN 

CALL  MUST1(R(I),C(I 
MR(I)-MBUY 
C(I)*C(I)+MR(I) 
S(I)«S(lj+MR(I) 
IC ( I ) -INT (  C  ( I )  + . 


C ( I ) , DLT ( I ) , SVL , MBUY ) 


MR(I)-MBUY 

C(I)*C(I)+MR(I) 

S(I)»S(lj+MR(I) 

IC(I)*INT(C(I)+.5) 

IS(I)*INT(S(I)+.5) 

EC(I)»EC1(I,II)+(MR(I)-DLT(I))*H(I) 


IIISMij 

C(I)-INT(C(I)  +  .! 


END  IF 

IF(KEY.EQ.2)THEN 

CALL  MUST2(R(I) , C(I) ,S(I) ,DLT(I) , SVL, MBUY) 

MR(I)*MBUY 

C(I)*C(I)+MR(I) 

S ( I ) *S ( I } +MR ( I ) 

IC(I)-INT{C(I)+.5) 

IS(lj*INT(S(I)+.5) 

EC(I )*EC1 (I , II )+(MR(I )-DLT(I) )*H(I) 

END  IF 

WRITE (7 , 505)1, R(I) ,MR(I) ,IC(I),IS(I),EC(I) 

SEC=SEC+EC(I) 

SHAVE* (TVC-SEC) /TVC 
WRITE ( 7 , 506 ) SEC , SHAVE 
STOP 

FORMAT ( IX, /, '**  SELECTED  OPTIONS  «'  /,1X,'(*  KEY  *',I2,2X, '*  SVL. 
&  =' ,F4.2,2X, '*  LEAD  TIME  TYPE  *' ,I2,2X, '*  DATA*' ,12, 

&2X , ' *  S.METHD*' ,12, IX, ')',//, 1  *  KEY*  1  :  SVL.  BY  PROB.  OF  NO  SHOR 
&TAGE  PER  REPLENISHMENT  CYCLE', /,'  2  :  SVL.  BY  FRACT.  OF  DEM 

A  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF  *  LTT*  1  :  CONSTANT 

&LEADTIME  2  s  RANDOM  (UNIFORM  W/  MEAN)  '  /,'  *  DATA=1  :  DATA 

&1  2  :  DATA2  3  »  DATA3  4  :  DATA4  *  SEARCH  METHOD  =1  i 

&  MODIFIED  SILVER  METHOD  2  s  HEURISTIC1, 

&///,20X, '***  INDEPENDENT  ORDER  POLICY  ***' 

&,//,lX, 'ITEM'  3X, 'LT  DEMAND ( LT . ) ' ,2X, ' EOO(N-INT. ) ' ,2X, ' (INTEGER) ' , 
&2X, 'REORDER' ,2X, 'EXP.  TVC(S) ' ,2X, 'S/O  COST’) 

FORMAT ( 13 , 2X, F6 . 2 , 2X, F4. 2) 

FORMAT ( 13 , 2X, F7 .3 , 3X, F5 . 2 , 2X, F5 . 2 ) 

FORMAT (//, 20X, ' ***  JOINT  ORDER  POLICY  ***',/ ,2X, ' ITEM'  5X, ' RHO' 
&,8X, 'REORDER' ,5X'CANBUY' ,5X. 'ORDERUP' ,5X, 'EXPECTED  COST') 

FORMAT ( IX, 13 ,4X, FT  . 3 , ' ( 1 , F5 . 3 , 1 ) ' , 3X , F8 . 3 , 6X , 14 , 7X, 13 , 4X, F9 . 3 , 3X, 
83) 

FORMAT (3 IX, 'SUM  OF  AVG.  ANNUAL  TVC($)=' ,F9.3) 

FORMAT  (  '  1  , / / , 15X , ' **  OPTIMAL  PAIR  (C,5)  ',/, 

&16X, 1  ITEM' ,7X, 'OPT.  C',7X ,'OPT.  S') 

FORMAT (16X, 13 ,9X,I4,9X, 14) 

FORMAT(2X,I3,4X,F6.4,7X,I4,9X,I4,7X,I4,7X,F10.4) 


FORMAT (45X, 'TOTAL  COST  ■  •  F10.4, ' (' ,F6.4, ' ) ' ) 

FORMAT ( 1 1 , / / , 3X , ' *  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT’ 
6, /,15X, ' (WHEN  LEADTIME  -  0) '  /.15X  '*  AT  ITERATION  »'  12,/, 

&1X  'ITEM'  4X  'RHO' ,4X, 'OPT.Cf ,2X, 'OPT. S' t3X, 'EXP. COST1) 

FORMAT ( 2X , 12 , 3X , F6 . 4 , 4X , 13 . 4X , I 3 , 4X . FS . 4 ) 

FORMAT (23X, 'TOTAL  •' ,F9.4, 1 ( 1 ,F6.4, ' ) ' ) 

END 

SUBROUTINE  POISON  (K£Y,DLT,Q,SVL,MX,M,MAX,CCDF1) 


PSO-l-SVL 
IF(KEY.EO.l)  THEN 

CALL  SeRVl ( DLT , PSO , MAX , M ) 

RETURN 
END  IF 

IF(REY.EQ.2)  THEN 
CALL  SERV2 (Q , DLT , PSO , MAX , M , CCDF1 ) 

RETURN 
END  IF 
END 

SUBROUTINE  MAXD(DLT,MAX) 

DIMENSION  PMF(750) 

1*1 

PMF ( 1 ) *EXP ( -DLT ) 

CDF=PMF ( 1 ) 

GO  TO  2 

PMF (I)* (DLT/ (I-l) )*PMF(I-1) 

CDF*CDF+PMF ( I ) 

IF(CDF.GE. 0.9995)  GO  TO  3 
1=1+1 
GO  TO  1 
MAX  =  1+1 
RETURN 
END 

SUBROUTINE  SERV1(DLT ,PSO ,MAX,M) 

DIMENSION  PMF (200),CDF( 200 ) , CCDF ( 200 ) 
PMF(1)=EXP(-DLT) 

CDF(1)*PMF(1) 

CCDF ( 1 )=l-CDr ( 1 ) 

CALL  MAXD(DLT,MAX) 

DO  1  1=2 ,MAX+1 

PMF(I)=(DLT/(I-1))*PMF(I-1) 

CDF ( I j=CDF ( I- 1 ) +PMF ( I ) 

CCDF(I)=1-CDF(I) 

DO  2  1=1 ,MAX+1 

IF(PSO.LT.CCDF(I) .AND. PSO. GT. CCDF (1+1) )GO  TO  3 
M=I 

RETURN 

END 

SUBROUTINE  SERV2 (Q , DLT , PSO , MAX , M , CCDF1 ) 

DIMENSION  PMF (200) , CDF (200) , CCDF (200) 

STON=Q*PSO 

PMF ( 1 )=EXP ( -DLT ) 

CDF(1)=PMF(1) 

CCDF ( 1)=1 .-CDF (1) 

CALL  MAXD(DLT,MAX) 

DO  1  1=2 ,MAX+1 

PMF(I)=(DLT/(I-1))*PMF(I-1) 

CDF(I)=CDF(I-1)+PMF(I) 

CCDF ( I )=1-CDF(I ) 

CONTINUE 

1=0 

IF(I . EQ.0)THEN 


c 

c 

c 


c 

c 


1 

3 

2 


C 

C 


1 

3 


PST-DLT 

ELSE 

PST*DLT-I-DLT*CDF(I)-*-I*CDF(I+l) 

END  XF 

IFgST.LT.STON)  THEN 

CCDF1“CCDF (M+l ) 

ELSE 
1=1+1 
GO  TO  2 
END  IF 
RETURN 
END 

SUBROUTINE  LEAD(TLT,TLT) 

CALL  LRND(719325,U, 1,1,0) 

TLYX2*TLT*U 

RETURN 

END 

SUBROUTINE  MUST1 (R, C,DLT , SVL ,MBUY) 

DIMENSION  PMF ( 400 ) , CDF ( 400 ) 

IC=INT(C+.5) 

RC=(1 . /R)**IC 
RHS=SVL/(R**IC) 

PMF ( 1 ) =EXP ( -DLT ) 

IF(PMF(l).LT.1.0E-70)  PMF(1)=0. 
CDF(1)=PMF(1) 

DO  1  1=2,400 

PMF(I)=(DLT/(I-1))*PMF(I-1) 

IF(PMF(I) .LT.1.0E-70)  PMF(I)=0. 
CDF(I)=CDF(I-1)+PMF(I) 

CONTINUE 

JS=-2 

SUM=0. 

KF=JS+2 
KL=JS+IC 
DO  2  IX=KF  KL 

SUM=SUM+PMF?IX+1 ) * ( ( 1 /R) **IX) 

TEST=RC*CDF ( KL+1 ) -PMF ( JS+2 ) - ( R** ( JS+1 ) ) *SUM 
IF(TEST.GE.RHS)  THEN 
MBUY=JS 
ELSE 
JS=JS+1 
GO  TO  3 
END  IF 
RETURN 
END 

SUBROUTINE  MUST2(R,C,S,DLT,SVL,MBUY) 

DIMENSION  PMF (400), CDF (400) 

IC=INT(C+. 5) 

RHS=(S-C+R*(1-R**C)/(1-R))*(1-SVL) 

PMF(lT=EXP(-DLT) 

IF(PMF(l).LT.1.0E-70)  PMF(1)=0. 

cdf(i)=pmf7i) 

DO  1  1=2,400 

PMF(I)=(DLT/(I-1))*PMF(I-1) 
IF(PMF(I).LT.1.0E-70)  PMF(I)=0. 

CDF ( I ) =CDF ( I - 1 ) +PMF ( I ) 

CONTINUE 
JS=-1 
SUM=0 . 

DO  2  IW=JS+1.JS+IC 
IF(IW.LT.l)  THEN 


I 

I 


c 

c 


VCDF-0 . 

eIF|  IW . EQ . 0 ) VCDF1*CDF ( 1 ) 
VCDF*CDF(IW) 


VCDF1*CD#(IW+1) 

END  IF 

SUM  =  SUM+(1/R**IW)*(DLT-IW-DLT*VCDF+IW*VCDF1) 

IF(JS.LT.l)  THEN 
VCDF*0 . 

IF(  JS . EQ . 0 ) VCDF1=CDF ( 1 ) 

ELSE 

VCDF=CDF(JS) 

VCDF1=CDF { JS+1 ) 

END  IF 

TESTaR**C*(DLT-JS-DLT*VCDF+JS*VCDFl) ♦ ( 1 -R) *R** ( C+ JS ) *SUM 
IF(TEST.LE.RHS)  THEN 
MBUY-JS 
ELSE 
JS=JS+1 
GO  TO  3 
END  IF 
RETURN 
END 

SUBROUTINE  OPTION  (ID,LT,K£Y,LTT,SVL,MSRCH) 

WRITE(6 , 50) 

READ(5,5l)  ID 
WRITERS, 53) 

READ(S,51)  KEY 
WRITE (6, 55) 

RE AD (5, 57)  SVL 
WRITe(6,52) 

READ(5,5l)  LT 
WRITE (6 , 54) 

READ? 5 ,51 )  LTT 
WRITE (6 ,56) 

READ (5 , 51 )  MSRCH 
RETURN 

* 

'50  FORMAT ( 1 1 '  '  *  ENTER  OPTION  FOR  DATA  FILE  (  1-  4)  I  1  « 

&DATA1  2  -  DATA2  3  =DATA3  4  =  DATA4 ' ) 

51  FORMAT(Il) 

52  FORMAT ( '  1 1  ,  '  *  ENTER  THE  LEADTIME  IN  MONTH  !  ' ,/,5X,  '1  =  1  MONTH 

&',/,5X,'3*  3  MONTH* ,/,5X, '6  *  6  MONTH* ) 

53  FORMAT ( '  1 1  '  *  ENTER  OPTION  FOR  THE  METHOD  OF  SVL.  MEASURE  ! 

&  1  *  PROB.  OF  NO  SHORTAGE  FER  REPLENISHMENT  CYCLE  ',/,4X,'  2 

&*  FRACTION  OF  DEMAND  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF  f) 

54  FORMAT  ( 1 1 '  ,  1  *  ENTER  THE  LEADTIME  TYPE  !  \/,5X,  *1  =  CONSTANT' 

&  1 ,/ , 5X, 1 2  *  RANDOM  (UNIFORM  W/  MEAN)') 

55  FORMAT ( '  1 '  ,  '  *  ENTER  THE  DESIRED  SERVICE  LEVEL  !  ' ,/,4X, ' (SHOULD 
&BE  IN  THE  FORM  OF  "F4.2"  j  EX.  ;  0.90  OR  0.99  ...ETC.)') 

56  FORMAT ( *  1 ' , *  *  ENTER  OPTION  FOR  SEARCH  METHOD  (1-2)!',/,'  1  = 

&MODIFIED  SILVER  METHOD  2  =  HEURISTIC') 

57  FORMAT (F4. 2) 

END 


1.  INPUT  FORMAT 


004 

1 

2 

3 

4 


050.00  0, 

290.000 
41.000 
77.000 
122.000 


20 

6.90 

1.20 

3.90 
2.30 


10.00 

10.00 

10.00 

10.00 


2.  SXEC  PROGRAM 


&TRACE  OFF 
&FN  «  PARA 

&FNO  ■  &CONCAT  OF  &FN  OUTPUT 

&TYPE  Do  you  need  to  compile  your  program  ?  (Y) 

&READ  VAR  &R  COMPILE 

&IF  &R  COMPlEE  NE  Y  &GOTO  -RUN 

-H  FORTVS  &FN 

&IF  &RC  EQ  0  &GOTO  -RUN 

&TYPE  Your  program  did  not  compile;  check  for  errors. 

&TYPE  Do  you  wish  to  view  the  program  LISTING  file?  (Y) 
&READ  VAR  &RSP1 

&IF  StRSPl  EQ  Y  BROWSE  &FN  LISTING  A 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y) 

&READ  VAR  &RESP1 

&IF  &RESP1  NE  Y  &EXIT  1 

&COMMAND  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP2 

&IF  &RESP2  EQ  Y  &GOTO  -H 

&EXIT  1 

-RUN 

FILEDEF  01  DISK  &FN  DATA1  A1 

FILEDEF  02  DISK&FN  DATA2A1 

FILEDEF  03  DISK  &FN  DATA3  A1 

FILEDEF  04  DISK  &FN  DATA4  A1 

FILEDEF  07  DISK  &FN  OUTPUT  A1  (LRECL  133 

LOAD  &FN  (START 

&IF  &RC  EQ  0  ScSKIP  9 

&TYPE  Your  program  did  not  run  correctly.-  check  for  errors. 
&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y) 

&READ  VAR  &RESP3 

&IF  &RESP3  NE  Y  &EXIT  2 

&COMMAND  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP4 

&IF  &RESP4  EQ  Y  &GOTO  -H 

&EXIT  2 

&TYPE  YOUR  OUTPUT  IS  IN  THE  FILE  &FN  OUTPUT  A 
&TYPE  Do  you  wish  to  BROWSE  your  output?  (Y) 

ScREAD  VAR  &RESP 

&IF  &RESP  EQ  Y  &C0MMAND  BROWSE  &FN  OUTPUT  A 
&TYPE  Print  your  output  file?  (Y) 

&READ  VAR  &RESP7 

&IF  &RESP7  EQ  Y  &C0MMAND  PRINT  &FN  OUTPUT  A 
-REDO 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y/N) 

&READ  VAR  &RE5P5 

&IF  &RESP5  EQ  Y  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP6 

&RESP56  =  &CONCAT  OF  &RESP5  &RESP6 
&IF  &RESP56  EQ  YY  &GOTO  -H 
&IF  &RESP6  EQ  Y  &GOTO  -RUN 
&EXIT 


3.  RUN  EXAMPLE 


Sara 

o  you  need-  to  compile  your  program  ?  (7) 

Is  FORTRAN  COMPILER  ENTERED.  11 >58 >02 

**MAIN**  END  OF  COMPILATION  1  ****** 

**POISON**  END  OF  COMPILATION  2  ****** 

**MAXD**  END  OF  COMPILATION  3  ****** 

**SERV1**  END  OF  COMPILATION  4  ****** 

**SERV2**  END  OF  COMPILATION  5  ****** 

**LEAD**  END  OF  COMPILATION  6  ****** 

**MUST1**  END  OF  COMPILATION  7  ****** 

**MUST2**  END  OF  COMPILATION  8  ****** 

**OPTION**  END  OF  COMPILATION  9  ****** 

VS  FORTRAN  COMPILER  EXITED.  11:58:10 


EXECUTION  BEGINS... 

*  ENTER  OPTION  FOR  DATA  FILE  (1-4)  ! 

1  a  DATA1  2  =  DATA2  3  =DATA3  4  =  DATA4 

1 

*  ENTER  THE  LEADTIME  IN  MONTH  ! 

1=1  MONTH 
3=3  MONTH 
6=6  MONTH 

1 

*  ENTER  OPTION  FOR  THE  METHOD  OF  SVL.  MEASURE  ! 

1  =  PROB.  OF  NO  SHORTAGE  PER  REPLENISHMENT  CYCLE 

2  =  FRACTION  0?  DEMAND  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF 

2 

*  ENTER  THE  LEADTIME  TYPE  ! 

1  =  CONSTANT 

2  =  RANDOM  (UNIFORM  W/  MEAN) 

1 

*  ENTER  THE  DESIRED  SERVICE  LEVEL  ! 

(SHOULD  BE  IN  THE  FORM  OF  "F4.2"  :  EX.  ?  0.90  OR  0.99  ...ETC.) 
0.99 

*  ENTER  OPTION  FOR  SEARCH  METHOD  (1-2)! 

1  =  MODIFIED  SILVER  METHOD  2  =  HEURISTIC 

1 

YOUR  OUTPUT  IS  IN  THE  FILE  PARA  OUTPUT  A 
Do  you  wish  to  BROWSE  your  output?  (Y) 

Irint  your  output  file?  (Y) 
n 

Do  you  wish  to  XEDIT  the  program  file?  (Y/N) 
n 

Do  you  wish  to  run  the  program  again?  (Y) 

R;  T=1.47/2.39  11:58:39 
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4.  OUTPUT  1  (  P2  -  0.99  ) 


**  SELECTED  OPTIONS  i 

(*  KEY  *  2  *  SVL.  -0.99  *  LEAD  TIME  TYPE  ■  1  *  DATA-  1  *  S.METKD-  1  ) 

*  KEY*  1  i  SVL.  BY  PROB.  OF  NO  SHORTAGE  PER  REPLENISHMENT  CYCLE 

2  «  SVL.  BY  FRACT.  OF  DEMAND  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF 

*  LTT-  1  «  CONSTANT  LEADTIME  2  :  RANDOM  (UNIFORM  W/  MEAN) 

*  DATA-1  s  DATA1  2  s  DATA2  3  :  DATA3  4  t  DATA4 

*  SEARCH  METHOD  *1  :  MODIFIED  SILVER  METHOD  2  :  HEURISTIC 


***  INDEPENDENT  ORDER  POLICY  *** 


ITEM 

LT  DEMAND (LT.] 

EOQ(N-INT. 

)  (INTEGER) 

REORDER 

EXP.  TVC($) 

S/O  COST 

1 

24.1671 

0.083 

158.800 

159 

25 

220.984 

1.983 

2 

3.4171 

0.083 

143.178 

143 

3 

34.383 

1.882 

3 

6.4171 

0.083 

108.840 

109 

7 

85.740 

3.499 

4 

10.1671 

0.083 

178.399 

178 

10 

82.217 

1.536 

SUM 

OF  AVG.  ANNUAL  TVC($) 

=  423.323 

*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  -  0) 

*  AT  ITERATION  :  1 


ITEM 

RHO 

OPT.C 

OPT.S 

EXP. COST 

1 

0.9942 

79 

146 

201.8498 

2 

0.9272 

31 

90 

21.7035 

3 

0.9650 

27 

35 

66.6963 

4 

0.9774 

45 

138 

63.6090 

TOTAL  *  353.3584(0.1641) 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  *  0) 

*  AT  ITERATION  :  2 


ITEM 

RHO 

OPT.C 

OPT.S 

EXP. COST 

1 

0.9972 

84 

152 

210.4767 

2 

0.9499 

36 

97 

23.4369 

3 

0.9768 

48 

89 

70.0377 

4 

0.9850 

71 

145 

66.8895 

TOTAL  *  370.3406(0.1240) 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  =  0) 

*  AT  ITERATION  :  3 

ITEM  RHO  OPT . C  OPT.S  EXP. COST 

1  0.9973  84  152  210.5555 

2  0.9464  36  96  23.1011 

3  0.9741  43  88  68.8164 

4  0.9835  70  143  65.8678 

TOTAL  -  368.3406(0.1299) 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  =  0) 

*  AT  ITERATION  i  4 

ITEM  RHO  OPT .C  OPT.S  EXP. COST 

1  0.9974  87  152  210.8846 

2  0.9469  36  96  23.1429 

3  0.9744  43  88  68.9584 

4  0.9835  70  143  65.9083 

TOTAL  =  363.3940(0.1286) 
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*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  *  0) 

*  AT  ITERATION  t  5 


ITEM 

RHO 

OPT.C 

OPT.S 

EXP. COST 

1 

0.9973 

87 

152 

210.8239 

2 

0.9467 

36 

96 

23.1271 

3 

0.9744 

43 

88 

68.9309 

4 

0.9835 

70 

143 

TOTAL  - 

65.8689 

368.75051 

ITEM 

1 

2 

3 

4 


RHO 

0.9973 

0.9467 

0.9744 

0.9835 


**  OPTIMAL 

PAIR  (C,S) 

** 

ITEM 

OPT.  C 

OPT.  S 

1 

87 

152 

2 

36 

96 

3 

43 

88 

4 

70 

143 

***  JOINT  ORDER 

POLICY  *** 

REORDER 

CANBUY 

ORDERUP 

EXPECTED  COST 

25 

112 

177 

211.9739 

-1  ' 

35 

95 

22.0671 

5 

48 

93 

67.8259 

7 

77 

150 

TOTAL  COST 

64.4122 
*  366.2788(0. 

5.  OUTPUT  2  (PI  *  0.95  ) 


**  SELECTED  OPTIONS  i 

(*  KEY  *  1  *  SVL.  -0.95  *  LEAD  TIME  TYPE  »  I  *  DATA-  1  *  S.METHD*  1  ) 

*  KEY-  1  i  SVL.  BY  PROB.  OF  NO  SHORTAGE  PER  REPLENISHMENT  CYCLE 

2  «  SVL.  BY  FRACT.  OF  DEMAND  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF 

*  LTT-  I  t  CONSTANT  LEADTIME  2  s  RANDOM  (UNIFORM  W/  MEAN) 

*  DATA-1  :  DATA1  2  s  DATA2  3  :  DATA3  4  :  DATA4 

*  SEARCH  METHOD  =1  :  MODIFIED  SILVER  METHOD  2  s  HEURISTIC 


***  INDEPENDENT  ORDER  POLICY  *** 


ITEM 

LT  DEMAND (LT.) 

EOQ(N-INT-) 

(INTEGER) 

REORDER 

EXP.  TVC($) 

S/O  COST 

1 

24.1671 

0.083 

158.800 

159 

33 

232.024 

15.113 

2 

3.4171 

0.083 

143.178 

143 

7 

35.343 

16.762 

3 

6.4171 

0.083 

108.840 

109 

11 

88.860 

22.051 

4 

10. 167 < 

0.083 

178.399 

178 

16 

84.977 

13.453 

SUM 

OF  AVG.  ANNUAL  TVC($) 

-  441.203 

*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  =  0) 

*  AT  ITERATION  s  1 

ITEM  RHO  OPT.C  OPT.S  EXP. COST 

1  0.9942  79  146  201.3498 

2  0.9272  31  90  21.7035 

3  0.9650  27  85  66.6963 

4  0.9774  45  138  63.6090 

TOTAL  -  353.8584(0.1980) 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  *  0) 

*  AT  ITERATION  :  2 

ITEM  RHO  OPT.C  OPT.S  EXP. COST 

1  0.9972  84  152  210.4767 

2  0.9499  36  97  23.4369 

3  0.9768  48  89  70.0377 

4  0.9850  71  145  66.8895 

TOTAL  =  370.8406(0.1595) 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  *  0) 

*  AT  ITERATION  :  3 

ITEM  RHO  OPT.C  OPT.S  EXP. COST 

1  0.9973  84  152  210.5555 

2  0.9464  36  96  23.1011 

3  0.9741  43  88  68.8164 

4  0.9835  70  143  65.8678 

TOTAL  =  368.3406(0.1651) 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  =0) 

*  AT  ITERATION  :  4 

ITEM  RHO  OPT.C  OPT.S  EXP. COST 

1  0.9974  87  152  210.8846 

2  0.9469  36  96  23.1429 

3  0.9744  43  38  68.9584 

4  0.9835  70  143  65.9083 

TOTAL  =  363.3940(0.1639) 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
{WHEN  LEADTIME  =  0) 

*  AT  ITERATION  i  5 

ITEM  RHO  OPT.C  OPT.S  EXP. COST 

1  0.9973  87  152  210.8239 

2  0.9467  36  96  23.1271 

3  0.9744  43  88  68.9309 

4  0.9835  70  143  65.8689 

TOTAL  *  368.7505(0.1642) 


**  OPTIMAL  PAIR  (C,S)  ** 

ITEM  OPT.  C  OPT.  S 

1  87  152 

2  36  96 

3  43  88 

4  70  143 


***  JOINT 

ORDER  POLICY 

ITEM 

RHO 

REORDER 

CANBUY 

ORDERUP 

EXPECTED  COST 

1 

0.9973 

32 

119 

184 

221.6339 

2 

0.9467 

4 

40 

100 

23.2671 

3 

0.9744 

9 

52 

97 

70.9459 

4 

0.9835 

13 

83 

156 

67.1722 

TOTAL  COST  =  383.0190(0.1319) 


APPENDIX  B 

*SIM':PROGRAM  TO  TEST  THE  MODEL 

VARIABLE  DEFINITIONS 

KEY  =  OPTIONS  FOR  ORDER  POLICY 

ISTA  =  OPTIONS  FOR  START  INVENTORY  LEVEL 

IPRINT=  OPTIONS  FOR  PRINTOUT 

NIT  =  NO. OF  ITEMS  IN  THE  GROUP  TO  BE  EVALUATED 

NS IM  =  SIMULATION  PERIOD 

SIMT  =  SIMULATION  INTERVAL  TIME  (1YEAR) 

DM  =  ANNUAL  DEMAND  RATE 

ME  =  ORDER-UP  POINT 

IEQ  =  ECONOMIC  ORDER  QUANTITY (INTEGER) 

MC  =  CAN  ORDER  POINT  (CAN  BUY) 

MR  =  REORDER  POINT (M1=IND. :M2= JOINT) 

TDA  =  TIME  OF  DEMAND  OCCURING  TIME 
TOA  =  ORDER  QTY.  ARRIVAL  TIME 
M  =  INDEX  OF  EVENT  FOR  DEMAND  OCCURED  ITEM 

IA  *  INDEX  OF  EVENT  FOR  OREDR  QTY.  ARRIVED  ITEM 

IET  =  INDEX  OF  EVENT  TYPE 

CT  =  MASTER  CLOCK  TIME  IN  YRS (THE  EARLIEST  EVENT  OCCURING  TIME) 

DA  =  TIME  INTERVAL  BETWEEN  2  SUCCESSIVE  DEMANDS 

ARV  =  UNIFORM  RANDOM  VARIABLE  FOR  LEADTIME  WITH  MEAN  "ULT" 

ULT  =  EXPECTED  LEADTIME 
IQ  =  ORDER  QTY. 

IT  =  INDEX  FOR  REPLENISH  TRIGGER  ITEM 

JT  =  INDEX  FOR  JOINT  REPLENISHED  ITEM 

I OH  *  ONHAND  INVENTORY 

I OHS  «  START  LEVEL 

UOH  =  UNIT  YEARS  ONHAND  QTY 

UBO  =  UNIT  YEARS  BACKORDER  QTY 

NBQ  =  NO.  OF  BACKORDERS  (QTY.) 

NTBQ  =  TOTAL  NUMBER  OF  BACKORDERS 

BOT  =  BACKORDER  TIME 

TBOT  =  TOTAL  BACKORDER  TIME 

NOD  =  TOATL  NUMBER  OF  ORDERS  FOR  EACH  ITEM 

NDM  =  TOTAL  NUMBER  OF  DEMANDS  FOR  EACH  ITEM 

NTDM  =  TOTAL  NUMBER  OF  DEMANDS  DURING  TMAX 

NJTO  =  NUMBER  OF  JOINT  REPLENISHMENTS 


UNIT  PRICE 

HOLDING  COST  RATE 

HOLDING  COST 

ORDER  &  REVIEW  COST 

INDIVIDUAL  ORDERING  COST 

JOINT  (GROUP)  ORDERING  COST 

BACKORDER  COST  PER  BACKORDER  UNIT 

BACKORDER  COST  (TIME  DEPENDENT  FACTOR) 

AVG.  ANNUAL  TOTAL  VARIABLE  COSTS 


DIMENSION  ISEED(IO) ,DM(10) ,MR(10) ,MC(10) .ME(10) ,IOH(10) ,ASVL(10) , 
&TDA( 100) , IT (10 , 100 )  . JT(10 , lOO) , IQ( 100),DA(100) ,T0A( 100 ) , IEQ( 10) , 
&H(10) ,OCI(10) ,UP(10) , SCU(IO) ,SCT(10) ,F(lO) .NJTO (100) , OCT (100) , 
&UOH(10,100) , AUOH( 100 ) ,  UBO( 10,100) , AUBO( 100) .NOD ( 10 , 100) ,KB(10,100 
&, ANOD( 100 ) ,NDM( 10,100) ,ANDM(100) ,NBQ(10,100) ,NTBQ(10, 100) ,KBT(10) 
ScBOT (10,100)  ,  TB0T(  10,100)  ,  FRC  (10,100),  FR0(  10,100)  ,  SVL  (10,100)  , 
&AFRO( 100 ) ,OC(10 , 100 ) , AOC(IGO) , COST ( 10 , 100 ) , TCOST(100 j ,NTOD( 100 ) , 
&HC0ST?10 , 100) , SCOSTi 10 , 100) , SCTU(10 ,100) ,SCTT( 10 ,100) , NTDM (100) , 
&AHCOST ( 100 ) , ASCTU ( 100 ) , ASCTT ( 100 ) , ABOT (10) , ABOQ ( 100 ) , ACOST ( 100 ) , 


^SCTU(lOO) ,ASCxT(100) ,ABOT(i0) ,ABOQ(100) ,ACOST(lOO 


non  non  non 


&THCOST ( 100 ) , TSCTU( 100  ) . TSCTT ( 100 ) ,CCTT ( 100 ) , Ml < 10 ) , M2 ( 10 ) . 

( 100 ) , TUBO ( 10(3 ) , TFRT ( 1 0 ) , TFRQ ( 1 00 ) , NTDMS ( 100 ) , NTODS ( 100 ) , 


&TU0H(1_ 

&ANB(10) 

COMMON  ICASE,ISEED,DM,DA,ARV,ULT,I,ISTP,MR,ME,KEY,IEQ 


C 

C 


DATA  INITIALIZATION 


SIMT 

TMAX 

IC 

ISD 

CT 

DT 

AT 

K 

N 

TVC 

NTJTO 

TANDM 

TANOD 

TA0C 

TCST 

STV 

ANJ 


*  1. 

=  0. 

*  0 

=  293715 
=  0. 

=  0. 

=  0. 

=  1 
*  0 
=  0. 

=  0 

*  0. 

=  0. 

=  0. 

=  0. 

=  0. 

=  0. 


-READ  PROCESS  OPTIONS  ----- 

CALL  OPTION (ID , RATE1 , RATE 2 , KEY , LTOPT ,NSIM , ISTA , IP ) 
-  READ  INPUT  DATA  - 


READ ( ID , 101 )NIT , OCG , ISV , MSV , LT 
ULT=LT/12. 

DO  401  1=1, NIT 


READ (ID , 100) DM (I ),M1(I),M2(I),MC(I),ME(I), IEQ(I ) ,OCI (I ) ,UP(I),F(I) 
IF(KEY.EQ.1)MR(I)=M1(I) 

IF (KEY.EQ. 2)MR(I )=M2(I) 


SCU(I )=RATE1*UP 


401 


I  1 

SCT(I ;  _ 

H£I)=UP(I)*F(I) 


RAT22*Up(l j 


CONTINUE 

WRITE(7, 198)ID, KEY, NSIM, ISTA, IP, LTOPT, MSV, ISV, ULT,LT 


DO  1  1*1, NIT 

ISEED(I)=ISD+I**3 
KBT(I)  =  0 
ANB(I)  =  0. 

TFRT ( I )  =  0. 
ABOT(I)  =  0. 
ASVL(I)  =  0. 
TOA(I)  =  10. **10 
TDA(I)  =  0. 

DO  2  1=1, NIT 
DO  2  J=1 , 100 
DA(J) 

iqU) 


UOH(I.J) 
TUOH(J) 
AUOH(J) 
COST ( I ,  J) 
TCOST( “* 
ACOST ( 
UB0(I 
TUBO( j; 
AUBO(J; 

nod7i , j) 

3D( 


:(j) 

■  j) 


NTOD(J) 


0. 

0 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0 

0 
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NTODS(J) 

ANCD(J) 

NDM(I.J) 

NTDM(J) 

NTDMS(J) 

ANDM(J) 

HBQ(i , J). 


* 


TFRQ(J! 
AFRO ( J 
OC(I, J 


oc  I,J 

OCT(J) 

OCTT(J) 

AOC(J) 


AOC(J)  =  0. 

NJTO(J)  =  0 

HCOST ( I , J)  =  0. 

THCOST(J)  =  0. 

AHCOST(J)  =  0. 

SCOST(I ,  J)  =  0. 

SCTU(I,J)  =  0. 

TSCTU(J)  =  0. 

ASCTU(J)  =0. 

SCTT(I ,  J)  =  0. 

TSCTT(J)  =  0. 

ASCTT(J)  =  0. 

KB  (I,J)  =  0 

IT  (I,J  =  0 

JT  (I.J)  =  0 

DO  3  1=1, NIT 

IF(ISTA.EQ.l . AND. KEY. EQ.l) THEN 
IOH(I)  =  IEQ(I) 

END  IF 

IF ( ISTA .EQ.l. AND . KEY . EQ . 2 )  THEN 
IOH(I)  =  ME ( I ) 

END  IF 

IF(ISTA.EQ.2)  THEN 
CALL  START 
IOH( I )=ISTP 
END  IF 
CONTINUE 


-  PRINTOUT  INPUT  DATA  - 


DO  4  1=1, NIT 

WRITE <7, 19?)  I,DM(I),IEQ(I),MR(I),MC(I),ME(I),IOH(I),ULT 
WRITE (7 ,97) 

DO  5  1=1, NIT 

WRITE(7 , 98)  I , UP ( I ) , F ( I ) , H ( I ) , OCG , OCI ( I ) , SCU ( I ) , SCT ( I ) 

IC  =  IC  +  1 

IF^IC.GE.  65. AND, IP. EQ.l. OR, IC.GE. 65. AND. IP. EQ. 2)  THEN 

WRITE (7, 200) 

END  IF 
CONTINUE  . 


-  START  SIMULATION  - 


IF  (IP .EQ.l .OR. IP. EQ. 2) THEN 
WRITE (7, 200) 

END  IF 


DO  777  ITER=1 ,NSIM 

IF  ( ITER . GT . 1 . AND . IP . EQ . 1 . OR . ITER . GT . 1 .AND . IP . EQ . 2 )  THEN 


W1' 


IC  *  0 

WRITE (7, 200) 

END  IF 

TMAX=TMAX+SIMT 

C 

20  IF(CT.GT.TMAX)  GO  TO  77 
N  *  N+l 
JOINT  =  0 
C 

C  -  DETERMINE  THE  INITIAL  DEMAND  OCCURANCE  TIMES  FOR  EACH  ITEM  - 

C 

IF(N.LE.l)  THEN 
DO  6  1=1, NIT 

ICASE=I 
CALL  DEMAND 
TDA(I)=DA(ICASE) 

6  U0H(  I , ITER)=U0H(I , ITER)+IOH(I)*DA(ICASE) 

C 

C  -  DETERMINE  THE  EARLIEST  EVENT  . 

C  (DEMMAND  OCCURED  /  ORDER  PLACED  /  ORDERED  QTY.  ARRIVED) 

C 

IF(TDA(1) .LE.TDA(2) )  THEN 
TD=TDA ( 1 ) 

ELSE 

TD=TDA(2) 

END  IF 

DO  7  1=1, NIT 
IF ( TDA ( I ) . LE . TD )  THEN 
TD=TDA(I) 

M=I 

END  IF 

7  CONTINUE 
C 

IET=0 
CT=TDA ( M ) 

I0H(M)=I0H(M)-1 

IF{MSV.EQ.1.AND.I0H(M) .EQ.-l)  KB(M,ITER)=KB(M,ITER)+1 
NDM ( M , ITER ) =NDM ( M , ITER ) + 1 
C 

C  -  SET  NEXT  DEMAND  OCCURING  TIME  FOR  THE  ITEM  : 'M'  - 

C 

ICASE=M 

CALL  DEMAND 

TDA(M)=CT+DA(ICASE) 

GO  TO  90 
END  IF 
C 

IF(TDA(1).LE.TDA(2))  THEN 
TD=TDA(1) 

ELSE 

TD=TDA ( 2 ) 

END  IF 

DO  8  1=1, NIT 

IF(TDA(I) .LE.TD)  THEN 
TD=TDA ( I ) 

M=I 
END  IF 

8  CONTINUE 
C 

IF(TOA(l).LE.TOA(2))  THEN 
TA=TOA( 1 ) 

ELSE 

TA=T0A ( 2 ) 

END  IF 

DO  9  1=1, NIT 

IF(TOA(I) .LE.TA)  THEN 
TA=T0A( I ) 

IA=I 
END  IF 
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9  CONTINUE 
C 

IF(TA.LE.TD)TKEN 
00  TO  93 
ELSE 
IBT»0 
CT*TD 
ICASE-H 
CALL  DEMAND  , 

TDA (M ) *CT+DA ( I CASE ) 

I0H(M)»I0H(H)-1 

IF(MSv.EQ.l .AND. IOH(M) .EQ.-l)  KB (M, ITER) *KB(M, ITER )+l 
NDM(M, ITER)aNDM(M, ITER)+i 
END  IF 
C 

C  -  DETERMINE  SHOULD  ORDER  - 

C 

90  IF(IOH(M).GT.MR(M))THEN 

UOH(M, ITER)*uOH(M, ITER)+IOH(M)*DA(M) 

IF(IP.£Q.1)THEN 

WRITE (7 ,501)CT, IET, M,UOH(M, ITER), (IOH(I) ,1*1,4) 

END  IF 

IF(IP.EQ.2)THEN 

WRITE(7 ,601 )CT, XET,M, UOH(M, ITER) , IOH(M) 

END  IF 

IC  *  IC  +  1 

IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  *  0 

WRITE (7, 200) 

END  IF 
GO  TO  20 
END  IF 


f(K.EO. 

IET*! 


1)THEN 


K*K+1 

IT(M,K)*M 

IF(IOH(M) .LT.O)  THEN 
NBQ(M,ITER)»-IOH(M) 
IFfKEY.EO.l)  IQ(M) 
IF(KEY.EQ.2)  IQ(M) 
IF (LTOPT.EQ . 1 )TOA(M 
IF(LT0PT.E0.2)THEN 
CALL  ARRIVE 
TOA(M)*CT+ARV 
END  IF 


KEY.EO.l)  IQ(M)  -  IEQ(M)+NBQ(M,ITER) 
KEY.EQ.2)  IQ(M)  ■  ME(M ) +NBO fn , ITER ) 
LTOPT.EQ. 1)T0A(M)  ■  CT  +  ULT 
LTOPT.EQ. 2) THEN 
LL  ARRIVE 


IF (NBQ  <M , I TER ) . EO . 1 ) TBOT (M , ITER ) »TBOT (M , ITER ) +BOT (M , ITER ) 
UBO (H , ITER) =UB0 ( ft , ITER ) +B0T ( M , ITER ) 

IF (IP. £0.1) THEN 

WRITE (7 , 502 ) CT , IET , M , NBQ (M , ITER) , BOT (M . ITER) , UBO (M , ITER ) , 
IQ(M) ,TOA(M) , (IOH(I) ,1*1,4) 

END  IF 

IF(IP.EQ.2)THEN 

WRITE (7, 602 )CT. IET, M,NBQ(M, ITER), BOT(M, ITER), UBO(M, ITER), 
IQ(M) ,TOA(M) , IOH(M) 

END  IF 

IC  *  IC  *  1 

IF  (XC.GE.65.AND.XP.EQ.1 .0R.XC.GE.65(AND.IP.EQ.2)  THEN 
IC  ■  0 

WRITE (7, 200) 

END  IF 
ELSE 

IF(KEY.EQ.l)  IQ(M)  ■  IEQ(M) 

IF(KEY.EQ.2)  IQ(M)«ME(MT-i6h(M) 

IF ( LTOPT . EQ . 1 ) TOA  <  M ) *CT+ULT 
IF (LTOPT.EQ. 2 )THEN 
CALL  ARRIVE 
TOA (M ) *CT+ARV 


vt  v-.*1  v  *  ■ 


non 


C 

C 


END  ZF 

IF<IP.EQ.1)THEN 

WRITET7 , 503 ) CT , IET ,M , IQ (M) , TOA(H) , ( IOH( I ) , I»1 , 4) 

END  IF 

XF(XP.EQ.2)THEN 

WRITE(7,o03)CT,IET,M,IQ(M) ,T0A(M) , IOH(M) 

END  IF 

IC  «  IC  +  1 

IF  ^  I C .  GE .  6  5 . AND .IP.EQ.1.0R.1C.GE.65. AND . IP . EQ . 2 )  THEN 

WRItI{7 ,200) 

END  IF 
END  IF 

NOD (M, ITER)  =  NOD (M, ITER )+l 
OC(M, ITER)  »  0C(M, ITER)+OCG+OCI (M) 


IF(KEY.EQ.l)GO  TO  20 
,2)GO  TO  92 


'.EQ.: 

IF(KEY.EQ.; 

END  IF 

-  DON'T  ORDER  UNTIL  PREVIOUS  ORDER  ARRIVES  — 

IF(M.EQ.IT(M,K) )THEN 
IET=0 
IF 
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CT 

IF (NBQ(M , ITER) . EQ . 1 )TBOT(M , ITER) =TBOT (M , ITER) +BOT (M , ITER) 

UBO (M, ITER )=UBO(M, ITER) +BOT(M, ITER) 

IF(IP.EQ.1)THEN 

WRITE  <7,504) CT  J)ET  ^  M^NBQ ( M , ITER ) , BOT ( M , ITER ) , UBO ( M , ITER ) , 

END  IF 

IF(IP.EQ.2)THEN 

WRITE (7 , 604 ) CT , IET , M , NBQ ( M , ITER ) , BOT ( M , ITER ) , UBO ( M , ITER ) , IOH ( M ) 
END  IF 

IC  *  IC  M 

IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  *  0 

WRITE (7, 200) 

END  IF 

ELSE 

UOH(M, ITER)=UOH(M, ITER)+IOH(M)*DA(M) 

IF(IP.EQ.1)THEN 

WRITE(7, 501)CT, IET, M,UOH(M, ITER), (IOH(I) ,1=1,4) 

END  IF 

IF(IP.EQ.2)THEN 

WRITE(7, 601 )CT, IET, M,UOH(M, ITER) ,IOH(M) 

END  IF 

IC  =  IC  +  1 

IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 

WRITE (7, 200) 

END  IF 
END  IF 
GO  TO  20 
ELSE 

DO  10  J=1,K 

IF(IQ(IT(M,J)).GT.O.OR.IQ(JT(M,J)).GT.O)GO  TO  92 
IET=1 
K=K+1 
IT(M,K)=H 

IF(IOH(M).LT.O)THEN 
NBQ (M, ITER)  =  -IOH(M) 

IFfKEY.EQ.l)IQ(M)  =  IEQ(M)+NBQ(M, ITER) 

IF  KEY.EQ.2)IQ(M)  =  HE (M)+NBQ(M, ITER) 

IF(LTOPT.EQ.llTOA(M)  =  CT  +  ULT 
IF(LTOPT.EQ.2)  THEN 
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non 


CALL  ARR 

toaM-c 

HD  IF 


RIVE 

CT+ARV 


IF(IP.E 

HRI 


BOT(M,ITER)  «  TOA(M)  -  CT 

IF (NBQ (M, ITER) . EO . 1 )TBOT?M ,ITER)=TBOT (M , ITER)+BOT (H , ITER) 
UBO(M^ITER^»UBO(fl, ITER)+BOT(H, ITER) 

WRITE ( 7 , 502 ) CT . IET.M, NBQ (M , ITER) , BOT (M. ITER) , UBO (M , ITER) , 
IQ(M),tOA(M7,(IOH(l),I=l,4} 


END  IF 
IF(IP.: 


(IP .EQ. 2) THEN 

WRITE (7 , 602) CT, IET ,M, NBQ (M, ITER) ,BOT(M, ITER )  ,  UBO  (M , ITER ) , 
IQ(M) ,TOA(M) , IOH(M) 


i  IQ(M) ,TOA(M) , IOH(M) 

END  IF 

IC  *  IC  +  1 

IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 

WRITE (7 ,200) 

END  IF 
ELSE 

IF(KEY.EQ.1)IQ(M)=IEQ(M) 

IF<KEY.EQ.2)IQ(M)=ME(M)-IOH(M) 

IF ( LTOPT . EQ . 1 ITOA ( M ) *CT+ULT 
I F ( LTOPT . EQ . 2 ) THEN 
CALL  ARRIVE 
TOA(M)*CT+ARV 
END  IF 

IF ( IP . EQ . 1 ) THEN 

WRITE(7, 503)CT, IET, M, IQ(M) ,TOA(M), (IOH(I) ,1=1,4) 

END  IF 

IF(IP.EQ.2)THEN 

WRITE (7 , 603 ) CT , IET , M , IQ (H) ,TOA(M) ,IOH(M) 

END  IF 

IC  =  IC  +  1 

IF  (IC.GE.65.AND.IP.EQ.1.CR.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 

WRITE (7, 200) 

END  IF 
END  IF 

NOD (M, ITER)  =  NOD(M,ITER)+l 
OC(M,ITER)  =  OC (M , ITER) +OCG+OCI (M) 

END  IF 

IF(KEY.EQ.l)GO  TO  20 


-  DETERMINE  THE  JOINTLY  REPLENISHABLE  ITEM  - 


91  DO  11  1=1, NIT 

IF  (IOH(I) .GT.MR(I) .AND.IOH(I) . LE.MC(I) .AND. IQ (I) .EQ.0)THEN 
JT(I.K)  =  I 
IQ(I)  =  ME(I)-IOH(I) 

IF ( LTOPT . EQ . 1 ) TOA  I ) =CT+ULT 
IF (LTOPT 


IF ( LTOPT . EQ . 1 ) TOA ( I ) =CT+ULT 
IF ( LTOPT. EQ. 2) THEN 
CALL  ARRIVE 
TOA ( I ) =CT+ARV 
END  IF 

NOD(I , ITER)=  N0D(I , ITER)+1 
OC(I , ITER)  =  0C(I , ITER)+OCI (I ) 
JOINT  =  1 


IF( IP .EQ. 1 .OR. IP . EQ . 2)THEN 
WRITE ( 7 , 506 ) I , IQ ( I ) , TOA ( I ) 
END  IF 


IC  =  IC  +  1 

IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
1C  =  0 

WRITE(7 ,200) 

END  IF 
END  IF 

11  CONTINUE 

NJTO(ITER)  *  NJTO(ITER)  + JOINT 


GO  TO  20 

92  IF<IOH(M).LT.O)  THEN 

NBQ{M,ITER)«-IOH(M) 

UBO(M,ITER)*UBO(H,ITER)+BOT(M,ITER) 

BOT (M, ITER) «TOA (M) -CT 

IF  IiPe'  1  tMn^  ‘ 1  >TB0T*M '  ITER)>TB0T<I('  ITER)+B°T(K, ITER) 
WRITE<7?5ol)CT,IET,M,NBQ(M,  ITER),  BOT(K,  ITER),  UBO(M,  ITER), 

&  (IOH(l) ,1*1,4) 

END  IF 

IF(IP.EQ.2)THEN 

WRITE(7,604)CT, IET, M,NBQ(M, ITER), BOT(M, ITER), UBO(M, ITER), IOH(M) 
END  IF 

IC  =  IC  +  1 

IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 

WRITE (7, 200) 

END  IF 

ELSE 

UOH(M,ITER)=UOH(M,ITER)+IOH(M)*DA(M) 

IF ( IP . EO . 1 ) THEN 

WRITEC7 ,501) CT, IET, M,U0H(M, ITER) , (IOH(I) ,1=1,4) 

END  iF 

IF(IP.EQ.2)THEN 

WRITE (7 , 601 )CT , IET,M, UOH(M, ITER) , IOH(M) 

END  IF 

IC  =  IC  +  1 

IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  *  0 

WRITE (7, 200) 

END  IF 
END  IF 
GO  TO  20 
C 

c  -  SET  NEW  ON  HAND  LEVEL  FOR  THE  ITEM  s'lA'  - 

C 

93  CT=*TA 
IET=2 

IOH(IA)=IOH(IA)+IQ(IA) 

DO  12  1=1, K 

IF7lA.EQ.IT(IA,I))THEN 

it7ia.i)=o 

IQ(IA)=0 

NTBQ  ( I A , ITER ) =NTBQ ( I A , ITER ) +NBQ ( I A , ITER ) 

BOT (I A, ITER) =0 . 

NBQ ( I A , ITER ) =0 
TOA(IA)*10>*10 
END  IF 

IF(IA.EQ. JT(IA, I) )THEN 
JT(IA,K)=0 
IQ ( I A ) =0 

NTBQ ( IA , ITER)=NTBQ ( IA , ITER ) +NBQ ( IA , ITER) 

BOTTiA,fTER)=0. 

NBQ (I A, ITER )=0 

toaIiaJ-io.^io 

END  IF 

12  CONTINUE 

IF(IP.EQ.1)THEN 

write ( 7, r - 

END  IF 

IF(IP.EQ.2)THEN 

WRITE ( 7, f - 

END  IF 

IC  =  IC  +  1 

IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 

WRITE (7, 200) 

END  IF 
GO  TO  20 


,  500)CT, IET , IA, (I0H(I) ,1=1,4) 
.  2)THEN 

,  600)CT, IET, IA, IOH( IA) 
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-  COMPOTE  INTER-SUMMARY  (ANNUAL)  - 

IF(IP.NE.4)THEN 
WRITE (7, 300) ITER 
END  IF 

DO  IS  1=1, NIT 


OCT (ITER ) =OCT ( ITER ) +0C ( I , ITER ) 
IF(1CV.EQ.1)THEN 

I F ( KB( I , ITER ) . GT . 0 ) THEN 

FRC(I,ITER)=KB(I,ITER)*1.0/NOD(I,ITER) 


ITT 

FRC(I,ITER)=0. 

END  IF 

SVL(I , ITER)=(1 . -FRC(I , ITER) )*100 . 

ELSE 

IF(NTBQ(I,ITER) .GT.O)THEN 

FRQ ( I , ITER ) =NTBQ ( I , ITER ) *1 . 0 /NDM ( I , ITER ) 

ITT 

FRQ(I,ITER)*0. 

END  IF 

SVL(I ,ITER)=(1 .-FRQ(I , ITER) )*100 . 

END  IF 
CONTINUE 
IF(IP.NE.4).THEN 
DO  16  1=1 , NIT 

WRITE (7 , 507)1 ,NDM(I , ITER) ,NOD(I , ITER) , UOH(I , ITER) ,UBO(I , ITER) , 
&NTBQ(I, ITER), TBOT(I, ITER), SVL(I, ITER), IOH(I) 

CONTINUE 


WRITE  (' 
WRITE  (' 
END  IF 


(7 . 508) NTDM(ITER) 

(7 . 509 ) NJTO( ITER) 


,NTOD(ITER) 


■  COMPUTE  COSTS  - 

DO  17  1=1, NIT 


SCOST ( I , ITER ) =SCTU ( I , ITER ) +SCTT ( I , ITER) 

COST (I , ITER ) =HCOST( I , ITER ) +SCOST ( I , ITER ) +0C ( I , ITER ) 
TCOST ( ITER ) =TCOST ( ITER ) +COST ( I , ITER ) 

IF(IP.NE.4) THEN 
WRITE(7 , 514) 

DO  18  1=1, NIT 

WRITE (7 ,515)  I , HCOST ( I . ITER) , SCTU ( I , ITER) ,SCTT(I , ITER) , 
&  OC(I , ITER) , COST \ I , ITER) 

WRITE (7, 516 )OCT( ITER ), TCOST (ITER) 

END  IF 

DO  19  1=1, NIT 

1  CONTINUE )=KBT<I)+KB<I'ITER> 

-  COMPUTE  FINAL-RESULTS  (MEAN  AND  STD.  DEVIATIONS)  - 

WRITE(7,510)NSIM 

DO  21  1=1, NIT 

DO  21  J*1 , ITER-1 

NTDMS ( I ) *NTDMS ( I ) +NDM ( I , J ) 

NTODS ( I ) =NT0DS ( I ) +NOD (I , J) 

TUOH( I )=TUOH ( I j  +UOH( I , J ) 

TUBO ( I ) *TUB0 ( I ) +UBO ( I , J ) 

TFRT ( I ) «TFRT ( I ) +TBOT ( I , J ) 

TFRQ ( I ) =TFRQ ( I ) +NTBQ ( I , J ) *1 . 0 
OCT? ( I ) =0CTT ( I ) +OC ( I , J ) 

THCOST ( I ) «THCOST ( I ) +HCOST (I , J) 

TSCTT ( I )«TSCTT ( I ) +SCTT ( I , J ) 


21 


22 

23 

24 

25 


26 


27 


97 


98 

100 

101 

198 


199 

200 


TSCTU< I )*TSCTU< I )+SC7U( I , J) 
DO  22  1*1, NIT 


ANDM 
ANOD 
AUOH 
AUBO 
ABOT 
ABOO(I 
ANB(I) 
AFRQ(I) 
IF(MSV.E 
IF (MSV .  E 


I 


■NfDMS(I)/TMAX 
■NTODS(I)/TMAX 
■TUOH7l)/TMAX 
bTUBO ( “ 

*TFRT( 

=TFRC 
*KBT( 


/TMAX 

,/TMAX 

/TMAX 


0/NTODS(I) 


)=TFRQ(I)/NTDMS(I) 

“*.1J  ASVL ( I )  =  ( 1 . -ANB ( I ) ) *100 . 


_  .2)  ASVL(l)*(l.-AFRQ(I))*100. 

AOC ( I ) =0£TT ( I ) / TMAX 
AHCOST ( I ) =THCOST ( I ) /TMAX 
ASCTU? I )=TSCTU( I ) /TMAX 
ASCTT ( I )=TSCTT ( I ) /TMAX 

ACOST ( I ) =AHCOST ( I ) +ASCTU ( I ) + ASCTT ( I ) +AOC ( I ) 

tvc=tvc+acost7i) 

DO  23  1*1, ITER- 1 
NT JTO=NT  JTO+N JTO ( I ) 

IF(NTJTO.GT.O)ANJ  =  NTJTO/TMAX 
DO  24  1*1, ITER-1 
STV=STV+ (TCOST  7 1 ) -TVC ) **2 
STDVC*SQRT ( STV/ ( ITER- 2 ) ) 

DO  25  1*1, NIT 

WRITE(7 ,511)1 ,ANDM(I) ,ANOD(I) ,AUOH(I) , AUBO(I) ,ABOQ(I) ,ABOT(I) , 

&  ASVL ( I ) 

DO  26  1*1, NIT 
TANDM=TANDM+ANDM ( I ) 

TANOD=TANOD+ANOD ( I ) 

TAOC*TAOC+AOC ( I ) 

CONTINUE 

WRITE ( 7 , 5 12 ) TANDM , TANOD 
WRITE(7,513)ANJ 
WRITE (7 ,517) 

DO  27  1*1, NIT 
WRITE (7, 515) 

WRITE?  / ,  §16)T, 

WRITE (7 , 518/STDVC 
STOP 

FORMAT (2X, ' 

&  PRICE'  2X, 'H.RA 
&/ , 36X , ' GROUP ' , 6X  ' IND . '  , 

FORMAT(2X,I2,4X,7(F6.2,3X) ) 

FORMAT (F7 .3 , 2X, 5(13 , 2X) ,3(F6.2,2X)) 

FORMAT* 13 ,3X, F6 .2 ,3X, I2,3X,2(I1,3X)) 

FORMAT ( / 7 , 5X , 1 **  INPUT  DATA  » '  / , 3X, / , 12X, ' *  USED  DATA  FILE  NO.  *' 
&, 13, /,12a,  '*  ORDER  POLICY  OPTION  = '  ,  13 , /  ,  12X,  ' *  STIMULATION  LIMIT 
&(YRS)  * 1 , 13  ,  / , 12X, 1  *  START  LEVEL  OPTION  = 1 , 13 , / , 12X, ' *  PRINTOUT 
&OPTION  *' ,13,/ ,12X, '*  LEAD  TIME  TYPE  = ' , 13 , / , 12X, /// ,7X, 1  *  ORDER 


.  I.AHCOST(I) , ASCTU(I) .ASCTT (I) ,AOC(I) ,ACOST(I) 
)TAOC,TVC 


, '  *  INPUT  DATA  FOR  COSTS* 1 ,//, '  ITEM' ,2X, 'UNIT 
'  ,2X, 'H.COST' ,3X  '  ORDER  COST' ,9X, 'B/O  COST' 

“  '  4X  'QTY. '  - 


lA,  '  UKUBH  » 

,5X, 'TIME’ ) 


&POLICY  OPTIONS' ,//,10X, '1  s  INDIVIDUAL  ORDER  POLICY' ,/, 10X, ' 2  :  JO 
&INT  ORDER  POLICY' ,//,7X, '*  START  LEVEL  OPTIONS ',//, 10X, ' 1  :  REGULA 

W/  E.O.Q  ;  JOINT  POLICY  W/  - - ‘ 

3INT  +  1; 


&R  (IND.  POLICY  , 

&2  t  RANDOM  7 REORDER  P03 

&OPTIONS' ,//,10X, '1  *  PRINT  ALL  ITEM' ,/, 10X, ' 2  :  PRINT  ONLY  EVENT 
&ITEM' , / , 10X, ' 3  s  PRINT  ONLY  SUMMARY' ,/,10X, ' 4  :  PRINT  ONLY  FINAL  R 
&ESULT' ,//,7X, '*  LEAD  TIME  TYPES  ' , // , 10X, ' 1  :  COSTANT  LEADTIME',/, 
&10X, ' 2  s  RANDOM  LEADTIME  (UNIFORM  W/  MEAN) ' ,// ,6X, ' *  SERVICE  LEVEL 


ORDER-UP  QTY .)',/, 10X , ' 
ORDER-UP  POINT)' ,//,7X, '*  PRINTOUT 
.  *  PR]" 


&  :  £' ,11, 

&I2 , '  MON. 


'*' ,13, IX  '%' ,5X, '*  AVG.  LEAD  TIME 


,  F6 . 3  '  YRS . ' , ' ( 


)',///, 3X , ' ITEM ' , 3X , ' YR  DEMAND'  4X, 'EOOf ,6X, ‘MUST  BUY' 
&3X , ' CAN  BUY' ,3X, 'ORD.  UP  QTY' , 3X ,' START  LEVEL 3X ,' EXP .  LT. ' ) 
. . ,4(I3,3X),3X,I3,7X,F7.3) 


FORMAT(3X, 12 ,6X,F7 .3. 5X 


FORMAT ( 

&  0  t 


**' ,//,20X, ' (EVENT  TYPE 
2  t  ORDER  QTY  ARRIVED)' 


'1'  //,35X,'**  SIMULATION  RESULTS 

DEMAND  OCCURED  1  *  ORDER  PLACED  _  _ 

&,//,5X, 'EVENT' ,4X, 'ITEM  #',2X,’BACK  ORDER 1 , 3X ,' UNIT  YfiS (QTY) ' , 4X , 
&ON  ORDER' ,10X, 'ON  HAND' ,9X, 'JOINT  REPLENISHMENT’ ,/ ,1X, 'TIME' ,2X, 

&'  ’"'TYPE  1  ,10X  'QTY.  '  ,2X,  'TIME'  ,3X,  'ON  HAND'  ,2X.  'B/O.  S'  ,3X,  'QTY.  '  ,2X, 
&' AR/T. ' j4X, ' 1  ,4X, ' 2 ' , 4X , '3' ,4X, '4' ,2X, ' ITEM#' ,2X, 'OD.  QTY. ' ,2X, '6 
&D.AR/T.  ) 
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no  n  on 


300 


500 
600 

501 
601 

502 

602 

503 

603 

504 

604 

506 

507 

508 

509 

510 


511 

512 

513 

514 


(ANNUAL)  ***'//. 

*  ITEM  # ' , 3X , 1 W  DEMANDS' ,.3X,f#  ORDERS' 


INTER  -  SUMMARY 


FORMAT? '1' ,///,20X,'***  INTI 

&• AT  THE  YEAR  « • ,2X,I3,//.3X,  **bll  TT  f  VA  ,  If  VbltfUlVii  t  §  TT  VfWbO 
&.10X,'  UNIT  YEAR  \7X, 'BACKOfoER' ,8X, 'SERVICE' ,6X,'F/  OH',/,40X, 

. T)',3X  '  B/0  ',2X,'0TY.',3X ,'TIME',5X, 'LEVEL(%)') 

'  f,Fd.4,2X,I2  3X  l2,48X,4?f4,lX)) 

'  ',F6.4  2X  12  3X,I2,54X,I4) 


'  , F6 . 4 , 2X , 12 , 3X , I 2 , 1 8X , F8 .2 , 22X , 4 ( 14 , IX) ) 
'  »?6‘4*2X,I2,3X,I2,1_8X,F8.2,28X.I4) 

• 3 / X IX § 


F6.2,3X,I3,1X, 


&f0N  RAND' 

FORMAT ( 

format! 

FORMAT ( ' 

FORMAT!' 

F0RMAT('  ' .F6.4;2X;i2;3X;i2;5X,l3,3X,F6. 

&  F6.3,1X,4?I4,1X)) 

FORMAT?'  ' , F6 . 4 , 2X , 12 , 3X , 12 , 5X , 13 , 3X , F6 . 3 , 11X , F6 . 2 , 3X , 13 , IX , 

&  F6.3,7X,I4) 

FORMAT!'  ' , F6 .4, 2X, 12 , 3X, 12 , 37X, 13 , IX, F6 .3 , 1X,4(I4, IX) ) 

FORMAT?'  ' , F6 .4, 2X, 12 , 3X, 12 ,37X, 13 , IX, F6 . 3 , 7X, 14) 

FORMAT?'  ' ,F6 .4,2X,I2,3X, 12 , 5X, 13 , 3X, F6 . 3 , 11X,F6 .2 , 14X,4(I4, IX) ) 
FORMAT?'  ' ,F6.3,2X,I2,3X,I2,5X,I3,3X,F6.3,11X,F6.2,20X,I4) 

FORMAT ?86X, 12 ,6X, 13 ,5X, F6 .3) 

FORMAT?5X,I2,9X,I4,/X,I3,9X,F9.4,2X,F6.3,2X,I3,3X,F6.3,5X,F6.2, 

&  9X, 14  5 

FORMAT (IX, 'TOTAL  = ' , 7X , 15 , 5X , 15) 

FORMAT? IX, 'TOTAL  NUMBER  OF  JOINT  REPLENISHMENT  =',I5) 

FORMAT? ' 1 ' .III ,20X, ' **  FINAL  RESULTS  (AVG.)  :' ,10X, 'DURING' ,13, 2X, 
& ' YEARS  ***  ' ,//,3X, 'ITEM  #' ,3X, '#  DEMANDS' ,3X, '#  ORDERS' ,9X, 'AVG 
&.  UNIT  YEAR' ,4X, 'AVG.  B/0' , 5X,' SERVICE'  ,/,40X, 

&' ON  HAND' , 3X, ' B/O.S ' ,4X, 'QTY. ' ,2X, 'TIME' ,3X, 'LEVEL' ,'(','%',')') 
FORMAT(5X, 12 ,6X, F6 .2 , 5X, F6 .2 , 10X, F7 .3 , IX, F7 .3 , IX, F6 .2 ,2X,F7 .3 , 
&2X,F6.2) 

FORMAT ( IX , 'TOTAL  AVG.=' , F7 .2 , 5X , F6 . 2) 

FORMAT? IX, ’AVG.  NUMBER  OF  JOINT  REPLENISHMENT  =',F6.3) 

FORMAT? IX, ///, 10X, '  *  COMPUTED  ANNUAL  COSTS  ($)',//, 2X, 

& ' ITEM ' ,3X, 'HOLDING' , 9X ,' BACKORDER ' ,9X, 'ORDER' ,6X, 'TOTAL' ,/, 22X, 'QT 
&Y. ' ,7X, 'TIME' ) 

FORMAT(2X, 12 , 5X, 5(F7 .3 ,4X) ) 

FORMAT ? IX, /,33X, 'TOTAL  *  ' , 2(F8.3 ,3X) ) 

FORMAT? IX, ///, 10X, '  *  COMPUTED  AVG.  ANNUAL  COSTS  ($)',//, 2X, 

&' ITEM' ,3X, 'HOLDING' ,9X, 'BACKORDER' ,9X, 'ORDER' ,6X, 'TOTAL' ,/,22X, 'QT 
&Y. ' ,7X, 'TIME' ) 

FORMAT ( IX, /,33X, 'STAND.  DEVIATION*  ' ,F8.3) 

END 

SUBROUTINE  DEMAND 

DIMENSION  ISEED(IO) ,DM(10) ,DA(100) ,MR(10) ,ME(10) ,IEO(10) 

COMMON  ICASE , ISEED , DM , DA , ARV , ULT , I , ISTP , MR , ME , KEY , IEQ 
CALL  LRND  (ISEED(ICASE) ,U, 1 , 1 ,0) 

DA<ICASE)=-(l./DM(ICASE))*ALOG(U) 

RETURN 


SUBROUTINE  ARRIVE 

DIMENSION  ISEED(IO) ,DM(10) , DA (100) ,MR(10) ,ME(10) ,IEQ(10) 
COMMON  ICASE , ISEED , DM , DA , ARV , ULT , I , ISTP , MR , ME , KEY , IEQ 
CALL  LRND(ISEED(I),U, 1,1,0) 

ARV*  2 . *ULT*U 

RETURN 

END 


SUBROUTINE  START 

DIMENSION  ISEED(IO) ,DM(10) ,DA(100) ,MR(10) ,ME(10) , IEQ(10) 
COMMON  I CASE , ISEED , DM , DA , ARV , ULT , I , ISTP , MR , ME , KEY , IEQ 
IF(KEY.EQ.l)  THEN 

CALL  LRND(ISEED(I) ,U, 1,1,0) 

IS*  (IEO(I)-(MR  I  +1  )*U 
ISTP*  IS+MR(I)+1 
RETURN 


c 

c 


30 

31 


END  IP 
!F(KEY.EQ.2 
CALL 


THEN 


ALL  £RND(ISEED(I) ,U, 1.1,0) 


ISTP* 
RETURN 
END  IF 
END 


D+l 


SUBROUTINE  OPTION (ID , RATE1 , RATE 2 , KEY , LTOPT ,NSIM , ISTA , IP ) 

WRITE (6 ,35) 

READ (5, 34)  ID 
WRITE (6, 39) 

READ (5 , 40 ) RATE 1 , RATE 2 
WRITE(6 ,31) 

RE AD (5, 34)  KEY 
WRITE(6 ,3/ ) 

READ (5, 34)  LTOPT 
WRITE(6 ,36) 

RE AD (5, 30)  NSIM 
WRITE(6 , 32) 

RE AD (5, 34)  ISTA 
WRITE(6,33) 

READ (5,34)  IP 
RETURN 
FORMAT (12) • 

FORMAT (  '  1 '  ,  '  *  ENTER  OPTION  FOR  ORDER  POLICY  (1  -  2)  !  ',/, 5X, 

&1  ■  INDIVIDUAL  ORDER  POLICY',/,  5X,'2  *  JOINT  ORDER  POLICY 


F  ') 

,/,5X, '1 


32  FORMAT ( 1 1 1 , 1  *  ENTER  OPTION  FOR  START  LEVEL  (1-2)  ! 

&=  REGULAR  ( IND .ORDER  POLICY  :  E.O.Q  /  JOINT  ORDER  POLICY  :  ORDER  U 
&P  QTY . ) ‘ , / , 5X , 1 2  *  RANDOM  BETWEEN  (MUST  ORDER  LEVEL+1 )  AND  (ORDER 
&UP  LEVEL)') 

33  FORMAT  (  1 1 '  ,  '  *  ENTER  OPTION  FOR  PRINTOUTS  -  4)  !  ',/,5X,'l  =  PRI 
&NT  ALL  ITEM' ,/,5X, '2  =  PRINT  ONLY  EVENT  ITEM' , / , 5X, ’ 3  *  PRINT  ONL 
&Y  SUMMARY  ',/,5X,'4  *  PRINT  ONLY  FINAL  RESULT') 

34  FORMAT(Il) 

35  FORMAT ( ' 1 ' , '  *  ENTER  OPTION  FOR  DATA  FILE  (1-4)  !  ',/, 

&5X , 1 1  *  DATA1 '  3X , 1 2  »  DATA2',3X,'3  *  DATA3',3X,'4  *  DATA4 1 ) 

36  FORMAT ( 1 1 ' , '  *  ENTER  SIMULATION  LIMIT  YEARS  ! 1 , / , 5X  ' 01  ■  5IMU 
&LATION  LIMIT  IS  1  YEAR'  / , 5X, ' 02  »  SIMULATION  LIMIT  IS  2  YEARS' 

&  ,/,7X  '  -  ' ,/,5X, '99  *  SIMULATION  LIMIT  IS  99  YEARS  ') 

37  FORMAT ( 1 1 ' , '  *  ENTER  THE  TYPE  OF  LEADTIME ( 1  -  2)  1', /,5X,'l  =  CON 
&STANT  LEADTIME ' , / , 5X , '2  =  RANDOM  (UNIFORM  WITH  MEAN)') 

39  FORMAT ( 1 1 ' , '  *  ENTER  RATE  OF  STOCKOUT  COST  TO  UNIT  PRICE  !',/,5X 
&,  '1.25  0.91  =  STOCKOUT  COST  PER  UNIT  IS  125%  OF  UNIT  PRICE 1 ,/, 17X, 
&' STOCKOUT  COST  PER  TIME  IS  91%  OF  UNIT  PRICE') 

40  F0RMAT(2(F4.2,1X) ) 

END 


1.  DATA  FORMAT  (  P2  =  0.99,  for  items  in  group  1  ) 


004  05C 

1.00 

99 

2 

1 

290.000 

025 

025 

112 

177 

159 

010.00 

006.90 

000.20 

041.000 

003 

-1 

035 

095 

143 

010.00 

001.20 

000.20 

077.000 

007 

005 

048 

093 

109 

010.00 

003.90 

000.20 

122.000 

010 

007 

077 

150 

178 

010.00 

002.30 

000.20 

2.  EXEC  PROGRAM 


&TRACE  OFF 
&FN  «  SIM 

&FNO  »  &CONCAT  OF  &FN  OUTPUT 

&TYPE  Do  you  need  to  compile  your  program  ?  (Y) 

&READ  VAR  &R  COMPILE 

&IF  &R COMPlEE  NE  Y  &GOTO  -RUN 

-H  FORTVS  &FN 

&IF  &RC  EQ  0  &GOTO  -RUN 

&TYPE  Your  program  did  not  compile;  check  for  errors. 

&TYPE  Do  you  wish  to  view  the  program  LISTING  file?  (Y) 
&READ  VAR  &RSP1 

&IF  &RSP1  EQ  Y  BROWSE  &FN  LISTING  A 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y) 

&READ  VAR  &RESP1 

&IF  &RESP1  NE  Y  &EXIT  1 

&COMMAND  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP2 

&IF  &RESP2  EQ  Y  &GOTO  -H 

&EXIT  1 

-RUN 

FILEDEF  01  DISK  &FN  DATA1  A1 

FILEDEF  02  DISK  &FN  DATA2  A1 

FILEDEF  03  DISK  &FN  DATA3  A1 

FILEDEF  04  DISK  &FN  DATA4  A1 

FILEDEF  07  DISK  &FN  OUTPUT  A1  (LRECL  133 

LOAD  &FN  (START 

&IF  &RC  EQ  0  &SKIP  9 

&TYPE  Your  program  did  not  run  correctly?  check  for  errors. 
&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y) 

&READ  VAR  &RESP3 

&IF  &RESP3  NE  Y  &EXIT  2 

&COMMAND  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP4 

&IF  &RESP4  EQ  Y  &GOTO  -H 

&EXIT  2 

&TYPE  YOUR  OUTPUT  IS  IN  THE  FILE  &FN  OUTPUT  A 
&TYPE  Do  you  wish  to  BROWSE  your  output?  (Y) 

&READ  VAR  &RESP 

&IF  &RESP  EQ  Y  &C0MMAND  BROWSE  &FN  OUTPUT  A 
&TYPE  Print  your  output  file?  (Y) 

&READ  VAR  &R&SP7 

&IF  &RESP7  EQ  Y  &COMMAND  PRINT  &FN  OUTPUT  A 
-REDO 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y/N) 

&READ  VAR  &RESP5 

&IF  &RESP5  EQ  Y  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP6 

&RESP56  =  &CONCAT  OF  &RESP5  &RESP6 
&IF  &RESP56  EQ  YY  &GOTO  -H 
&IF  &RESP6  EQ  Y  &GOT0  -RUN 
&EXIT 


63 


3.  RUN  EXAMPLE 


sin 

Do  you  need  to  compile  your  program  ?  (Y) 

Is  FORTRAN  COMPILER  ENTERED.  12 :11s 11 

**MAIN**  END  OF  COMPILATION  1  ****** 

**DEMAND**  END  OF  COMPILATION  2  ****** 

**ARRIVE**  END  OF  COMPILATION  3  ****** 

**START**  END  OF  COMPILATION  4  ****** 

**OPTION**  END  OF  COMPILATION  5  ****** 

VS  FORTRAN  COMPILER  EXITED.  12:11:23 


EXECUTION  BEGINS . . . 

*  ENTER  OPTION  FOR  DATA  FILE  (1-4)  ! 

1  =  DATA1  2  =  DATA2  3  *  DATA3  4  =  DATA4 

1 

*  ENTER  RATE  OF  STOCKOUT  COST  TO  UNIT  PRICE  !  ' 

1.25  0.91  a  STOCKOUT  COST  PER  UNIT  IS  125%  OF  UNIT  PRICE 
STOCKOUT  COST  PER  TIME  IS  91%  OF  UNIT  PRICE 

0.40  0.01 

*  ENTER  OPTION  FOR  ORDER  POLICY  (1-2)  ! 

1  =  INDIVIDUAL  ORDER  POLICY 

2  =  JOINT  ORDER  POLICY 

1 

*  ENTER  THE  TYPE  OF  LEADTIME (1  -  2)  ! 

1  »  CONSTANT  LEADTIME 

2  =  RANDOM  (UNIFORM  WITH  MEAN) 


*  ENTER  SIMULATION  LIMIT  YEARS  ! 

01  =  SIMULATION  LIMIT  IS  1  YEAR 
02  =  SIMULATION  LIMIT  IS  2  YEARS 

99  =* SIMULATION  LIMIT  IS  99  YEARS 
50 

*  ENTER  OPTION  FOR  START  LEVEL  (1-2)  ! 

1  =  REGULAR (IND. ORDER  POLICY:  E.O.Q  /  JOINT  ORDER  POLICY:  OD.UP  QTY. ) 

2  *  RANDOM  BETWEEN  (MUST  ORDER  LEVEL+1)  AND  (ORDER  UP  LEVEL) 


*  ENTER  OPTION  FOR  PRINTOUT (1  -  4)  ! 

1  =  PRINT  ALL  ITEM 

2  =  PRINT  ONLY  EVENT  ITEM 

3  »  PRINT  ONLY  SUMMARY 

4  =  PRINT  ONLY  FINAL  RESULT 
4 

YOUR  OUTPUT  IS  IN  THE  FILE  SIM  OUTPUT  A 
Do  you  wish  to  BROWSE  your  output?  (Y) 

frlat  your  output  file?  (Y) 
n 

Do  you  wish  to  XEDIT  the  program  file?  (Y/N) 
n 

Do  you  wish  to  run  the  program  again?  (Y) 

R;  T=4.63/5.82  12:11:54 


4.  OUTPUT  1  (individual  order  policy) 


**  INPUT  DATA  i 

*  USED  DATA  FILE  NO.  ■  1 

*  ORDER  POLICY  OPTION  »  1 

*  SIMIMULATION  LIMIT (YRS)  *  50 

*  START  LEVEL  OPTION  =  2 

*  PRINTOUT  OPTION  =  4 

*  LEAD  TIME  TYPE  =  1 


*  ORDER  POLICY  OPTIONS 

1  s  INDIVIDUAL  ORDER  POLICY 

2  :  JOINT  ORDER  POLICY 

*  START  LEVEL  OPTIONS 

1  ;  REGULAR  (IND.  POLICY  -  E.O.Q  ;  JOINT  POLICY  -  OD.UP  OTY.) 

2  :  RANDOM  (REORDER  POINT  +  1;  ORDER-UP  POINT) 

*  PRINTOUT  OPTIONS 

1  :  PRINT  ALL  ITEM 

2  s  PRINT  ONLY  EVENT  ITEM 

3  :  PRINT  ONLY  SUMMARY 

4  :  PRINT  ONLY  FINAL  RESULT 


*  LEAD  TIME  TYPES 

1  :  COSTANT  LEADTIME 

2  :  RANDOM  LEADTIME  (UNIFORM  W/  MEAN) 


*  : 

SERVICE  LEVEL 

s  P2=  99  % 

*  AVG.  LEAD  TIME:  0.083  YRS. 

(  1  MON.) 

ITEM 

YR  DEMAND 

EOO 

MUST  BUY  CAN  BUY  ORD.  UP 

START 

EXP.  LT 

1 

290.000 

159 

25 

112  177 

65 

0.083 

2 

41.000 

143 

3 

35  95 

45 

0.083 

3 

77.000 

109 

7 

48  93 

38 

0.083 

4 

122.000 

178 

10 

77  150 

60 

0.083 

*  INPUT 

DATA  FOR  COSTS; 

ITEM 

UNIT  PRICE 

H.RATE 

H.COST  ORDER  COST 

B/O  COST 

GROUP  IND .  Q 

TY. 

TIME 

1 

6.90 

0.20 

1.38 

50.00  10.00  2 

.76 

0.07 

2 

1.20 

0.20 

0.24 

50.00  10.00  0 

.48 

0.01 

3 

3.90 

0.20 

0.78 

50.00  10.00  1 

.56 

0.04 

4 

2.30 

0.20 

0.46 

50.00  10.00  0 

.92 

0.02 

**  FINAL  RESULTS  (AVG, 

.)  s  DURING  50 

YEARS 

*** 

ITEM 

# 

#  DEMANDS 

#  ORDERS 

AVG.  UNIT  YEAR  AVG. 

B/O 

SERVICE 

ON  HAND  B/O.S  OTY. 

TIME 

LEVEL (% 

1 

292.24 

1. 

84 

79.151  0.025  2792 

0.009 

99.00 

2 

41.58 

0. 

30 

69.217  0.008  0.34 

0.006 

99.18 

3 

77.00 

0. 

72 

54.064  0.004  0.46 

0.003 

99.40 

4 

120.50 

0. 

68 

88.760  0.014  0.84 

0.006 

99.30 

TOT  AVG. 

531.32 

3. 

54 

AVG. 

NO. 

OF  JOINT 

REP;  0. 

000 

iMllMH 


*  COMPUTED  AVG.  ANNUAL  COSTS  ($) 


HOLDING 

BACKORDER 

ORDER 

TOTAL 

QTY. 

TIME 

109.228 

8.059 

0.036 

110.400 

227.723 

16.612 

0.163 

0.002 

18.000 

34.777 

42.170 

0.718 

0.003 

43.200 

86.090 

40.830 

0.773 

0.007 

40.800 

82.409 

TOTAL 

=  212.400 

430.999 

STAND. 

DEVIATION  = 

51.942 

5.  OUTPUT  2  (joint  order  policy) 


**  INPUT  DATA  i 

*  USED  DATA  FILE  NO.  »  1 

*  ORDER  POLICY  OPTION  *  2 

*  SIMIHULATION  LIMIT (YRS)  =  50 

*  START  LEVEL  OPTION  =  2 

*  PRINTOUT  OPTION  =  4 

*  LEAD  TIME  TYPE  *  1 


*  ORDER  POLICY  OPTIONS 

1  :  INDIVIDUAL  ORDER  POLICY 

2  s  JOINT  ORDER  POLICY 

*  START  LEVEL  OPTIONS 

1  :  REGULAR  (IND.  POLICY  -  E.O.Q  ;  JOINT  POLICY  -  OD.  UP  QTY.) 

2  s  RANDOM  (REORDER  POINT  +  1;  ORDER-UP  POINT) 

*  PRINTOUT  OPTIONS 


1  s  PRINT  ALL  ITEM 

2  :  PRINT  ONLY  EVENT  ITEM 

3  :  PRINT  ONLY  SUMMARY 

4  ;  PRINT  ONLY  FINAL  RESULT 

*  LEAD  TIME  TYPES 


1  :  COSTANT  LEADTIME 

2  :  RANDOM  LEADTIME  (UNIFORM  W/  MEAN) 


*  SERVICE  LEVEL:  P2=  99  %  *  AVG.  LEAD  TIME:  0.083  YRS.(  1  MON.) 


ITEM 

YR  DEMAND 

EOO 

MUST  BUY 

1 

290. 00C 

159 

25 

2 

41.000 

143 

-1 

3 

77.000 

109 

5 

4 

122.000 

178 

7 

CAN  BUY 

ORD.  UP 

START 

EXP.  LT 

112 

177 

71 

0.083 

35 

95 

28 

0.083 

48 

93 

32 

0.083 

77 

150 

50 

0.083 

*  INPUT  DATA  FOR  COSTS: 


ITEM 

UNIT  PRICE 

H.RATE  H.COST  ORDER  COST 

B/O  COST 

GROUP 

IND, 

QTY. 

TIME 

1 

6.90 

0.20  1.38 

50.00 

10.00 

2.76 

0.07 

2 

1.20 

0.20  0.24 

50.00 

10.00 

0.48 

0.01 

3 

3.90 

0.20  0.78 

50.00 

10.00 

1.56 

0.04 

4 

2.30 

0.20  0.46 

50.00 

10.00 

0.92 

0.02 

**  FINAL  RESULTS  (AVG, 

.)  : 

DURING  50  YEARS 

*** 

ITEM 

# 

#  DEMANDS 

#  ORDERS 

AVG.  UNIT  YEAR 

AVG.  B/O 

SERVICE 

ON  HAND 

B/O.S 

QTY.  TIME 

LEVEL  (% 

1 

292.24 

1.94 

75.843 

0.021 

2.56  0.007 

99.12 

2 

41.58 

0.60 

56.035 

0.000 

0.00  0.000 

100.00 

3 

77.00 

1.18 

51.541 

0.001 

0.04  0.001 

99.95 

4 

120.50 

1.06 

30.055 

0.003 

0.30  0.001 

99.75 

TOT  AVG. 

;  531.32 

4.78 

AVG. 

NO. 

OF  JOINT 

REP;  1.620 

I>UMH 


*  COMPUTED  AVG.  ANNUAL  COSTS  <$) 


HOLDING 

BACKORDER 

ORDER 

TOTAL 

104.663 

OTY. 

TIME 

7.066 

0.030 

113.400 

225.159 

13.448 

0.000 

0.000 

6.000 

19.448 

40.202 

0.062 

0.001 

12.800 

53.065 

36.825 

0.276 

0.001 

13.600 

50.703 

TOTAL 

*  145.800 

348.375 

STAND. 

DEVIATION  = 

21.055 
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